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Abstract 

We will establish the connection between the Lorentz covariant and so- 
called single-time formulation for the quark Wigner operator. To this 
end we will discuss the initial value problem for the Wigner operator 
of a field theory and give a discussion of the gauge-covariant formula- 
tion for the Wigner operator including some new results concerning the 
chiral limit. We discuss the gradient or semi-classical expansion and 
the color and spinor decomposition of the equations of motion for the 
Wigner operator. The single-time formulation will be derived from the 
covariant formulation by taking energy moments of the equations for 
the Wigner operator. For external fields we prove that only the low- 
est energy moments of the quark Wigner operator contain dynamical 
information. 
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I. INTRODUCTION 



Ultra-relativistic heavy ion experiments nowadays reach extreme beam energies 
of up to 200 GeV/A at CERN. The compressed nuclear matter becomes so hot and 
dense that a Quark-Gluon-Plasma (QGP) is expected to form [|T|. The quest for this 
new state of matter |^ motivates the following work. It will deal with the description 
of ultra-relativistic heavy ion collisions in the framework of a transport theory for the 
underlying degrees of freedom, namely quarks and gluons. So far Quantum Chro- 
modynamics (QCD) as the gauge theory for quarks and gluons has been successful 
in describing equilibrium properties (including linear response functions to external 
perturbations) through (i) perturbation theory in the framework of thermal field the- 
ory (TFD) [^3, especially in connection with the Hard Thermal Loop resummation 
scheme of Braaten and Pisarski , and (ii) (non-perturbatively) through numerical 
Monte Carlo simulations of lattice gauge theory Since, however, the formation 
of a QGP is, at least in the early dynamical stages, governed by non-equilibrium 
processes we also need a kinetic theory which can deal with non-equilibrium phase 
space dynamics. For this the proper framework is relativistic transport theory. 

Classical relativistic kinetic theory has been used in practice for relativistic 
heavy ion collisions in the form of numerical cascade codes, including mean field as 
well as quantum statistical scattering effects for fermions (Pauli suppression) ||10| , p 



and bosons (stimulated scattering) |l2| , pr3| . A classical relativistic transport theory 
for colored degrees of freedom was developed in [0,|T^. The basis for quantum 
transport theory is the Wigner formalism |T6|-pO|]; through a systematic gradient 



expansion it can be reduced to classical kinetic theory 



Quantum transport equations for gauge theories require a gauge covariant defi- 
nition of the Wigner operator [ p!^P^p2| , p3[ . Furthermore, proper inclusion of field 



theoretic retardation effects is facilitated by a Lorentz covariant formulation [IIU 



On the other hand, if one wants to solve the kinetic equations as an initial value 
problem in time, a Lorentz covariant formulation leads to conceptual difficulties. 
These will be shortly reviewed in Sec. 0. This led Bialynicki-Birula, Gornicki, and 
Rafelski (BGR) ||2^ to propose a so-called single-time formulation for relativistic 
quantum transport theory which breaks explicitly the Lorentz covariance of the 



theory but can be solved as an initial value problem |2J-|26|. 

Our goal in this paper is to establish the connection between the Lorentz covari- 
ant and single-time formulations of quantum transport theory. The present work 



extends the previous study by Zhuang and Heinz to the case of non-Abelian 



gauge theories. We concentrate on the dynamics of fermion fields in the background 
of a (classical or fluctuating) gauge field. We start from the Lorentz covariant formu- 
lation of relativistic quantum transport theory, discussing its equations of motion 
and their semiclassical expansion as well as their color and spinor decomposition 



(Sec. p^ . Most of this reviews (sometimes in more elegant notation) previous work 
which is relevant for our analysis here. We do, however, point out that quite gener- 
ally transport and mass shell equations need to be solved only for 8 of the 16 spinor 
components, either for the scalar, pseudoscalar, and tensor densities or for the vec- 
tor and axial vector densities. The remaining 8 components can then be obtained 
from constraints (Sec. fTTT]). This was apparently not realized before. Also the 



color decomposition in Sec. [1IID| has not been presented before. Some new results 



are also contained in Sec. [Ill G| on the chiral (massless) limit and in Sec. [111H| on 



the classical limit. In particular we show that in the chiral limit the kinetic equa- 
tions for the right-handed and left-handed components decouple which considerably 
simplifies the structure of the covariant transport equations. 



In Sec. 13 we derive the single-time form of the transport equations by tak- 



ing energy moments of the covariant equations. This method was established in 



Ref. for simpler theories and is here extended to the case of general (abelian 
or non-abelian) gauge theories. The covariant transport equations are equivalent 
to an infinite hierarchy of coupled single-time transport equations for the energy 
moments of the covariant Wigner operator. We construct this hierarchy and discuss 
its practical truncation for applications. The general moment hierarchy is written 
down for the Wigner operators and their equations of motion which, after taking 
ensemble expectation values, would generate the whole coupled BBGKY hierarchy 
of kinetic equations for 1-body, 2-body, etc. phase-space distribution functions. 

In the present paper the truncation of the energy moment hierarchy is discussed 
explicitly only for the simple limit of external gauge fields (Hartree limit). This 
corresponds to Vlasov-type transport equations without collision terms. In this limit 
we find a simple and exact truncation scheme ( "exact" meaning that no additional 
approximations like e.g. a semi-classical gradient expansion are neccessary) in which 
only one dynamical single-time transport equation for the lowest energy moment 
must be solved while all higher order moments can be obtained from the solution 
via constraints. This explains the success of the BGR approach (who only 
studied the lowest energy moment and never looked at the full moment hierarchy) 



for external field problems and at the same time provides a basis for generalizing 
it to include back-reaction and collision effects. In the last section of the main 
text, Sec. [iV D| , we write these equations out explicitly for QED to leading order 
of a semiclassical gradient expansion to recover the well-known Vlasov equations 
supplemented by spin precession terms. 

A short summary of the achievements and an outlook on the remaining open 
problems is given in Sec. 0. The Appendix contains helpful formulae and other 
technical details. 

II. THE INITIAL VALUE PROBLEM 

The Wigner operator for a field theory with second quantized fields \E', \E'^ is 
given by the Fourier transformation of the density matrix 

^(a;+,x_) = ^(a;+)*"^(x_), x± = x±y/2 (2.1) 

over the relative coordinate ?/ = — x_ 

W{x,p)^ /^e-^-^^(x + |,x-|). (2.2) 

X = (x+ +x_)/2 is the center of mass coordinate, and p is the canonically conjugate 
momentum to y. 

Using the translation operators exp{—y ■ dx/2) (acting to the right) and exp{y ■ 
9|,/2) (acting to the left) to rewrite the density matrix, 

^ (x + |, a; - f ) = * (a; + f ) ^^(x - f ) = ^(x)e^ (^^-^-/2)^^(x) 

= <i^{x)e'^-y^\x), n='-{dx-dl), (2.3) 



the connection with the classical phase space density becomes apparent: 



= ^{x)6{p-tt)^^{x). (2.4) 



TT = i{dx — dl)/2 is the momentum operator. The Wigner operator (more exactly: 
its ensemble expectation value) thus is a quantum mechanical and off-shell general- 
ization of the classical phase space density 

fix, p, t) = E - W) ^(P - P^ii))^ (2-5) 

i 

which gives the probability density to have a particle with (on-shell) momentum p 
at space-time coordinate {x,t) (where x = (x+ + x_)/2). In contrast to the classical 
phase-space density /, W{x,p) = {W{x,p)) and its energy average W{x,p,t) = 
J dpoW{x,p) are in general not positive definite, but can become negative at small 
phase-space distances. Only if averaged over sufficiently large phase-space volumes 
AV = A^x A^p ^ (27r)'^ with some positive definite weight function g{x,p), 

f{x,p,t)= / ——fw{x-x\p-p\t)-g{x\p'), (2.6) 

the resulting f{x, p, t) is positive definite and can be interpreted as a classical phase- 



space distribution [16 . 



A. Initial value problem and energy moments 

For a covariant quantum field theory the density operator ) and, in 

general, all 2-point functions contain two time parameters, ti = xo + yo/2 and 
t2 = — yo/2. The occurrence of two time parameters is necessary for a proper 



description of signal propagation and retardation effects. It poses, however, concep- 
tual problems for a transport theoretic approach because the corresponding covariant 
Wigner density cannot be initialized with the help of asymptotic field configurations 



at t = — oo [^. To see this let us rewrite the covariant Wigner function (see (p.2|) ) 
as follows: 

(2.7) 



Obviously, due to the Fourier transformation over relative times yo, the covariant 
Wigner function at any fixed time Xq requires knowledge of the fields \E', \E'"'" at 
all times. at — >■ — c>o can thus not be initialized by giving the fields only at 
t = — oo. 

The equal-time (or single-time) Wigner function 

P^^Iwf ("^ + 2 ' ^o) (a; - f , xo) ) , (2.8) 

which uses the fields at equal times and thus depends only on a single time parameter, 
has clearly no such problem. Initializing W3 at xo = —00 by vacuum solutions for 
the fields, it has been used in P3 to calculate the pair production rate in QED in 



a constant external electric field, reproducing the correct Schwinger rate On 
the other hand, by comparing ( |2.7] ) and ( |2.8D one sees that the equal-time Wigner 



function is simply the energy integral of the covariant Wigner function W4: 

/oo 
dpoW,{x,p,po). (2.9) 
-00 

As such it contains less dynamical information than W4 since the complete off-shell 
structure of PV4 is averaged over. 



While the formulation of the transport theory as an initial value problem thus 
seems to require a single-time language, it can clearly not be based on W^, alone. 
Retardation and memory effects in the covariant approach are reflected in the off- 
shell behaviour of the covariant Wigner function W4 which should not be thrown 
away by simply averaging it over pq. Instead, one must study, in addition to Ws 
(which is the zeroth po-^ioment of W4) also the higher po-nioments of the covariant 
Wigner function: 



(27r)3 

[9S)v^(a; + f,a:o)] [^^-'^^^^ - f, ^o)] ) . (2.10) 



fc=0 

X 



The full information on is equivalent to knowing all (and not just the lowest)of 
its po moments, and initializing them at t = —00 requires the knowledge of all time 
derivatives of the fields \1', \E'^ at t = —00 (or, equivalently, to knowing the fields for 
all times). 

This is clearly impractical. It follows from these considerations that a formula- 
tion of covariant transport theory as an initial value problem for equal-time energy 
moments of the covariant Wigner function must be accompanied by a truncation pre- 
scription which allows to restrict the description to a finite number of such energy 
moments. 

The strategy of the present paper will thus be to first derive from the underlying 
field theoretic Lagrangian a set of covariant equations of motion for the covariant 
Wigner operator W4 and then, by taking po-moments of these covariant equations, 
set up the coupled hierarchy of equal-time equations for the energy moments of W4 
2^j3y]. We will then study the truncation of this hierarchy. 



In general, truncating the moment hierarchy will entail nontrivial approxima- 
tions. However, if the fermion fields evolve under the infiuence of an external gauge 
field, i.e. they do not interact with each other, neither directly via scattering nor 
indirectly via backreaction of their electric current on the external field, their time 
evolution is uniquely fixed once their values at t = — oo have been specified. (The 
Dirac equation is a first order differential equation in time.) It should thus be 
sufficient to specify initial conditions for W3 only, but not for any higher energy mo- 
ments. This can only be true if the time evolution of the system in the equal-time 
framework can be described completely in terms of W3 only. This turns out to be 
correct in the sense that, in the external field limit, only for a genuine equal-time 
transport equation must be solved while all higher energy moments of W4 can be 
obtained from the solution of this equation via constraints. This will be discussed 



in more detail in Sec. IV. 



B. Formulation with orthogonal polynomials 



In Ref. 1^31 suggested to construct the equal-time moments of the covariant 

Wigner density W4 from orthogonal polynomials of po rather than from the powers of 
Pq used in Eq. ( p.lO[ ). Orthogonal polynomials have the advantage that the complete 



function W4 is easily reconstructed from the set of moments using the orthogonality 
relations. Their disadvantage is that, if they are defined over the full Pq interval 
(—00, 00), they require additional weight functions which fall off sufficiently rapidly 
at ±00. Only if one restricts po to a finite interval [—a, a] one can use orthogonal 
polynomials without weight function, the Legendre polynomials. This was done 



in Ref. |^. The price for doing so was the technical requirement that the Wigner 
density IV4 falls off as a function of Pq sufficiently fast near ±00 that its contributions 
outside a large but finite interval can be neglected. 

If one wants to avoid this technicality and rather work with orthogonal polyno- 
mials on the infinite po interval, one must use additional weight functions g{po) in 
order to obtain suitable orthogonality relations. We will now show that this does not 
really lead to a proper single-time formulation. Using the reduced energy variable 
z = po/E with E = \^p'^ + m'^, the general moment expansion then reads 



00 -j^ 

W{x,p,Po) = J2—^k[X,p)Uk[Z) 



Wk[X,p) 



00 
00 



dz g{z) Uk{z) W{x, p, Ez) 
dzg{z) Uk{z) ui{z) . 



(2.11) 
(2.12) 
(2.13) 



Inserting the definition ( p.2| ) of the Wigner operator we get 



Wk{X,p) 



dzg{z) 
d^y 



d^y 



Uk{idyjE)e-'P-y 



Uk{idyjE) / dzg{z)e 



-lEzyo 



dz g{z)e 



-iEzyo 



Uki-idyjE)g{x+,x^) , (2.14) 



with partial integration in the last line, assuming that g{x^, x_) vanishes identically 

at |xo+|, |xo_| 



■ 00. However, Eq. (|2.14|) only then strictly corresponds to a 
formulation in terms of only a single time parameter if 



dz g{z) e 



-iEzyo 



6{Ez)=6{yo). 



(2.15) 



Since for orthogonal polynomials on the full po interval g{z)^l this is not possible. 



a strict single-time formulation cannot be achieved in this way. 



If we insist on a true single-time formulation there are thus only two possibilities: 
(1) One tries to keep the computational advantages of orthogonal polynomials. In 
this case one must restrict po to a finite interval and assume that outside that interval 
the Wigner function vanishes [^. (2) One does not want to place such a drastic 



restriction on W4 but allows for, say, exponential tails in the Wigner function at 
large po- In that case one should use the moments ( p.lOj ) with simple powers of po 
which are not orthogonal on each other. We will here take the second choice. 

Let us remark that the moment expansion of the covariant Wigner operator needs 



not be restricted to energy moments. In [31| a systematic expansion of the Wigner 



operator into moments of all four components of the momentum p was given. One 
can show that such an expansion corresponds to the semiclassical h expansion or, 
equivalently, to a gradient expansion around x. We will explain this in more detail 



in Sec. ^ when we discuss the connection between the expansion in energy moments 



and a temporal gradient expansion around xq. 

C. Green functions and the Wigner operator 

Related to the initial value problem of non-equilibrium dynamics in field theory 
is the doubling of the Hilbert space between zero and finite temperature field theory 

The connection between the non-equilibrium Wigner formalism and the two- 
point Green functions of a field theory in thermodynamic equilibrium is given via 
the density operator |l33| : 



iG{x+,x_) = (0|T^(x+,a;_)|0). (2.16) 



T is the path ordering operator for a given time contour C. The most perfect formal 
agreement []TRpD|,p2| between vacuum field theory, finite temperature equilibrium 



thermodynamics and non-equilibrium dynamics is obtained if one chooses as the time 
contour C the "Closed Time Path" (CTP) [|3 which goes from t = — cxDtot = -|-oo 
infinitesimally above the real time axis, then returns to t = —oo infinitesimally 
below the real time axis, and (for equilibrium thermodynamics) finally proceeds 
vertically to — oo — i/5 (where j3 = 1/T) in order to ensure the KMS condition, i.e. 
the periodicity of the thermal Green functions in imaginary time. The last piece 
is missing in non-equilibrium situations in which case the Green functions do not 
satisfy the KMS condition. 

In the CTP formalism the 2-point Green functions can be split into four pieces 
according to the four possibilities to distribute the two time arguments on the upper 
and lower parts of the real-time contour. In matrix notation one writes |T8|,20,32 



G{x,y) 



(2.17) 



\G>ix,y) G%x,y) ] 
The covariant Wigner function W4 is the 4-dimensional Wigner transform (Fourier 
transform with respect to y) of G'^{x + y/2,x ~y/2). 

In thermodynamic equilibrium formal identities, in particular the KMS condi- 
tion, allow to express all four components of the 2x2 matrix ( p.l7| ) in terms of a 



single real spectral density In non-equilibrium situations, the absence of the 
KMS condition implies that at least two components of ( |2.17| ) are independent and 
probably should be kept for a complete dynamical description. While G'^ and 
contain also information on the virtual excitations of the vacuum, and G^ (which 



are again related by the spectral density) describe exclusively the dynamics of the 
real excitations in the medium [^U|,^. Here we will be mostly interested in the case 
of the external fields where knowledge of the dynamics of resp. is sufficient. 



III. LORENTZ-COVARIANT FORMULATION 



In the following we will discuss the gauge theories QED and QCD |^2|,^. We 



will consider only the fermion Wigner operator and treat the photon or gluon degrees 
of freedom through their classical gauge fields. Problems connected with a suitable 



gauge choice will not be discussed here, the reader is referred to |22,36,37| 



A. Gauge-covariant Wigner operator 

For a gauge covariant theory the Wigner operator must transform covariantly 
under gauge transformations. This is achieved through a gauge covariant definition 



of the density operator p!^j22 



q(x + 1,x-1) = ^(x)e^-^'(")/2 ® e-y'^-^^^^^ix), (3.1) 

with the covariant derivative 

D^{x) = d^- tgA^ix) (3.2) 

instead of the partial derivative in the translation operator of Eq. ( p^.3| ). The direct 
product in Eq. (|3T| ) is over spinor and (in the case of QCD) color indices. Both 
the fields ip'^, a = 1,2, ... , Nc, and the covariant translation operators are group 
elements of SU (Nc) albeit in different representations. 
With ( p.l|) the Wigner operator reads 



(3.3) 



Similar manipulations as in Eq. (U) show that using; the covariant instead of the 
partial derivative results in p being the physical or kinetic momentum rather than 



the canonical one |14,22|. The gauge covariant translation operators in (3J,0) can 



be rewritten as link operators ||2^ 



+ =^(x + f)f/(a; + f,x) ®f/(x,a;-f,x)V'(a;-f) (3.4) 



with 



(3.5) 



where the integral in the exponent is taken along the straight path connecting b and 
a, and P denotes path ordering of the operators. The link operators are unitary and 



fulfill 



U\a,b) = U-\a,b) = U{b,a) 



U{a,b)U{b,a) = U{a,a) = 1. 



(3.6) 



The density operator can then be rewritten as 



(x + |, a; - l) = ^ (x + l) f/ (x + |, x) ® f/ (x, X - l) (x - l) . (3.7) 



In the following we will use the conventions from ||2^ for an operator O operating 
on the color and spinor indices of g: 



6 ^ (x + f , X - f ) =7/^ (x + f ) [/ (x + f , x) ® dU (x, X - f ) (x - f ^ 
^ (x + f , X - f ) 0=^ (x + f ) [/ (x + f , x) O ® t/ (x, X - f ) ^/^ (x - f ) . (3.^ 



Color octet objects in QCD will be denoted as matrices, O = Oata, with ta being the 
generators in the fundamental representation of the gauge group. They are given in 
Appendix |^ together with useful color trace formulae. Occasionally we will also use 
their adjoint representation Ta with 

{Ta)bc = -i'hfabc (3.9) 

and norm 

tr(T„Tfe) = ?,h^5ab. (3.10) 

With our normalization the ta have eigenvalues ~ | proportional to the physical 
color-spin of quarks, while the Ta have eigenvalues ~ h corresponding to the color 
spin of gluons. 

^From the Wigner operator one obtains the operators for the fermionic baryon 
(charge) and color currents via 

jl-y°\x) = -tr [^,i,{x)^{x)) = I d^p tr (l3x37.W^(a;,P)) , (3.11) 
jlTix) = -tr (^^,^^{x)4j{x)j = J d'p tr (^^^^W{x,p)j . (3.12) 

The baryon current projects onto the color singlet, the color current onto the color 
octet contributions of the Wigner operator. In analogy we will later on perform a 
color decomposition for the equations of motion of the Wigner operator. 
As spinor matrices g and W transform according to 

^t(x + f,x-f) =70^(x-f,x + f)7° (3.13) 

with exchanged arguments x± and 

W\x,p) = 7°#(x,p)7°. (3.14) 



B. Equations of motion for the Wigner operator 

The basis for the dynamics of strong interactions is the Lagrange density of QCD 
(since flavor quantum numbers don't matter for our purpose we consider only one 
quark flavor): 

£qcd = ihc^{x)-f''D^ix)ij{x) - mc^lfij - -l-F^,{x)F^'- {x). (3.15) 

ip{x), il){x) = ip'^'j^ are the spinor fields, 

D^{x) = d^- YcM^)^ (3-16) 

is the covariant derivative (see Eq. (|3.2|) ), m is the bare mass, and F^^{x) is the 
field strength tensor. The latter is defined through the commutator of the covariant 
derivative, 

F^,{x) = -^[D^,D,] = {d,A,{x)) - {d,A,{x)) - '£[A^ix),A,{x)], (3.17) 

and thus satisfies the Jacobi identity 

[Daix), F^,{x)] + [D,{x), Fa^{x)] + [D^{x), F^^x)] = 0. (3.18) 

Since for the semiclassical and gradient expansions we will later need to count powers 
of h we displayed them here once explicitly. Until further notice we will, however, 
from now on return again to natural units with h = c = 1. 
Inserting the Dirac equation and its adjoint, 

(i-f^'D^ix) - m)^(x) = 0, (3.19) 
i^{x){i-f''Dl{x) - m) = 0, (3.20) 



(with the hermitean adjoint of the partial derivative defined as before as acting 
to the left) into the definition ( ^.71 ) of the density operator and, following P^p5 



pulling the derivatives out in front by using the derivative rules for link operators as 
given in Appendix 0, one finds the following equations of motion for the covariant 



fermion Wigner operator 



2mW{x,py^^'[[tD^{x) +Pt^-9J ds (1 - 2s) + isdp)id';)w{x,p) 

+W{x,p)[-iD^{x) + 3 ds (1 + 2s) ["lF,^(a; + isdp)id^ 



(3.21 

2mW{x,py[(-iD^{x)+p^ + g y°^^ ds (1 + 2s) ["^^.^(x + isdj,)id;)w{x,p) 
+W{x,p)(iD^{x)+p^ - g ^ds{l- 2s) ^'^F,^{x + tsdp)id;) 



(3.22) 

Here, for notational convenience, momentum derivatives standing to the right of the 
Wigner operator are defined in the sense of partial integration as 

w{x,p)d;' ...d;'^ = {-ifd;*^ . ..d;'W{x,p) . (3.23) 

The operator 

[^]F,^(z(s)) = U{x, z{s))F,^{z{s))U{z{s), x) with z{s) = x + sy (3.24) 

is called Schwinger string and connects the field strength tensor at point z gauge 
covariantly with the point x along a straight line path. The Schwinger string is the 
essential quantity to describe the gluonic degrees of freedom as non-local operators. 
To get a more compact notation we define generalized non-local momentum and 



derivative operators 11^ and A^, respectively^ (here we display h and c explicitly for 
later use): 



= + ^ /°^^ dsishdl ^""^F^^ix + ishdp) , 

C J-l/2 



The equations of motion for the Wigner operator are then 

2mW{x,p) = Y ({H;,, W{x,p)} + t[A^, W{x,p)] 
2mWix,p) = {{U^,W{x,p)}-t[A^,W{x,p)]) 7^ 

With Hi = n^, (^A^)t = tA^ and [n^,7°] = [A^,7°] = we find 

7°(n^±zA,J V = U^±tA^, 



(3.25) 
(3.26) 



(3.27) 
(3.28) 



(3.29) 



which means that Eqs. (|3.271) and ( |3.28|) are adjoint to each other: 7° (|3.27|) 7° 
( p.28|) "^. Adding and subtracting the two equations gives 



4mW{x,p) = {7^ {n^, W{x,p)}} + z[7^ [A^, W{x,p)]] , 
= [7^ {n^, W{x,p)}] + ^{7^ [A^, Wix,p)]} . 



(3.30) 
(3.31) 



This corresponds to a separation into real and imaginary parts. The generalized mo- 
mentum operator is always working as an anti-commutator, the generalized deriva- 
tive operator always as a commutator on W. The only difference between the 



*Our notation here differs somewhat from that of Ref. |35]; it is optimized for a simuL 



taneous description of non-abelian and abelian gauge interactions, while in only the 
abelian case of QED with classical external fields was studied. In that limit A and 11 as 
defined here reduce to the definitions given in II, see Eqs. ( ^ , |3^ ) below. 



two equations is that for the spinor structure commutators are replaced by anti- 
commutators. Defining therefore 

= (3.32) 

we can combine the two equations as 

[2ml, W{x,p)]± = [Y, {n^, W{x,p)}]± + 2[7^ [A^, W{x,p)]U. (3.33) 

These are the equations of motion for the Wigner operator for particles with spin 
1/2. 

The equations are closed by the Yang-Mills-equations for the field strength tensor. 
These couple to the Wigner operator for the fermion fields via the current ji,: 

[Df^ix), F^,{x)] = j^ix) = tail taluW{x,p). (3.34) 

The trace sums over spinor and color degrees of freedom and contains an integral 
over momentum space. 

Eqs. ( p. 33 ) and ( 3.34 ) contain the full dynamical information of the system on 



the operator level. 

C. Gradient expansion 

The classical limit of these equations is obtained via a gradient expansion of the 
non-local operators 11 and A in Eqs. ( p.25| , pT^BD . It is generated by a Taylor expan- 



sion for the Schwinger string around x in powers of ishdp. As such it is automatically 
also an expansion in powers of h; since, however, the color decomposition discussed 
in the following subsection introduces additional powers of h, we will refer to the 



semiclassical expansion in the sense of Wigner and Kirkwood only as the "gradient 
expansion" . As we will see it is actually an expansion in covariant derivatives. It 
thus preserves gauge covariance. 



The Taylor expansion of the Schwinger string is given by 



k=0 



y=0 



(3.35) 



The connection between the expansion of the Schwinger string into local operators 
[{dp ■ dyY ^'^^F^nix + y)\y=Q and a gradient expansion in covariant derivatives D^{x) 
can be shown from the expansion of the Schwinger string in an exponential as derived 



in 22 : 



+ ishdp) = exjp[ishdp'Daix)]Ff,i,{x) 

CO -| ^ 



k=0 



(3.36) 
(3.37) 



T>a{x) is defined as the covariant derivative in the adjoint representation |p2| : 



V^{x)F^,{x) = [D^{x),F^,{x)]. 



(3.38) 



Comparison between Eqs. ( p.35|) and ( |3.37|) gives the following identity (see also 
Appendix for a direct calculation from of the l.h.s. of Eq. ( p.38| ) for n = 1): 



{dp-dyT t"^F,^(a; + y) 



y=0 



dp-D{x)r,F,,{x) 



(3.39) 



Here we defined a generalized commutator through 



[A", B] = B, 



(3.40) 



or explicitly 



E 

fc=0 



n' 



(3.41) 



n times 

Using Eqs. ( |3.35y3.39|) in Eqs. ( |3.25| , p.2^ ) and performing the s-integration we find 



for the non-local operators 



a. = ^^.-^|E(- 



k=0 



2 J {k + 



l(d,-D{x)r,F^,{x)]d;, (3.42) 



ih\ k + 1 



2c 



k=0 



2 J {k + 2)\ 



l{d,-D{x))\F^,{x)jd;. 



(3.43) 



Using the definition ( p.l7|) for the field strength tensor and in ( p.43|) the identity 



= [dp ■ D{x),Pf,] 



(3.44) 



we finally get 



= D^{x) 



oo 1 

Eyli-i9p-Dix))',nD,ix) 



oo 1 _ L 

U, = p, + 2j:^l{-'idp-D{x)r, 



k=l 



(3.45) 
(3.46) 



Eqs. ( p.45| , |37^ ) identify and 11^ as non-local generalizations of the covariant 
derivative hD^^ and the momentum p^. 

The (hermitean) operator for the expansion about the classical limit is thus given 

by 



A = 9f ■ iD^ix) 



(3.47) 



with h as the expansion parameter. The superscript W means that the momentum 
derivative acts on the momentum dependence of the Wigner operator. The super- 
script F on the covariant derivative reminds us that it acts only on the field strength 
operator following it, see Eqs. ( p.42| , |3l43|) . From Eqs. ( p.45| , p746D or (|3.44|) we find 
that A^ is one order higher in A and h than 11^. 



Since, at least for QED, the measurable classical fields are the field strengths 
F^u{x) we will use for the gradient expansion Eqs. ( |3.42| , ^.4^ ) rather than 
Eqs. (^^g^. 

The validity of the gradient expansion requires 



h/\ = hd^ ■ iD^{x) < 1, 



(3.48) 



i.e. that the Wigner function is sufficiently smooth in momentum space and the 
field strengths vary sufficiently slowly (in a covariant sense) in coordinate space. 
The corresponding length scales on which these functions vary must satisfy 



(3.49) 



Inserting the expansions ( p.42| , p^43D together with the Ansatz 

oo 

for the Wigner operator into the equations of motion ( p.33| ) we obtain 



(3.50) 



[2ml, 



(3.51) 



n-l 



r 



9_ 

2c,t^J/c + l)! 



i 



± 



-A)^F.,(x)la;,^y("-'=-l 



With Table | for the zeroth and first order terms in the non-local operators we get 



[2ml,W^(o)]±= [7^2p^#(o)]± 



(3.52) 



to zeroth order and 



7^2p,w^(i)-^{F.,(x)a;,iyw} 

9 



D,(x)-^F.,(x)9;,#w 



(3.53) 



to first order in h. 

We will examine the resulting classical limit more closely in the context of the 
spinor decomposition for the Wigner operator since there we will find an easy con- 
nection to the classical mass shell condition and transport equation. 



D. Color decomposition 



For QCD it is useful to decompose the Wigner operator into its (observable) 



color singlet and its (unobservable) color octet contributions |]14 



W{x,p) 
Ws{x,p) 



W,{x,p)l + W''{x,p) t" 



-ti W{x,p), 



2 



^tit" W{x,p) 



(3.54) 
(3.55) 
(3.56) 



The color decomposition of the Schwinger string |3^ can be calculated to any fixed 
order of the gradient expansion from the explicit form 



Wf,.(x + zshd,) = Y: ^-^liT^d^a^ + ^d;Al{x)tT , Ft^{x)t% (3.57) 

n=0 ^ '^^ 



Combining the two decompositions we find to zeroth order in the gradient expansion 



singlet: [ml - 7%, W^i°^]± = 0, 
octet: [ml - 7>^, H^(°)"]± = 0, 



(3.58) 
(3.59) 



and to first order 



singlet: [ml - 7>^, W^^^]^ + ^[7^ d^Wi% 



24c V2 



;[7^5;[f;^(x),#(°)i]± + [i^.d;{F:^{x),w^'^^}y 



-^[7^[i;:(x),iv(°)i]^) (3.60) 

octet: [ml - 7%, W^^^'']± + ^[7^, S^iy^^^"]^ = 



8c ' V2 

-'£[rAA',{x),W^''>^}Uy (3.61) 

Here we used the trace formulae given in Appendix ^ and the relation [AB, C] = 
A[B,C] + [A,C]B. 

These equations exemplify a general result which to prove we will leave to the 
reader: At any given order n of the gradient expansion, the "leading" color singlet 
and octet contributions W^"'\W^"'^ decouple from each other; the color singlet and 
octet sectors couple only via the lower order components W^"'~^\ 

Note that the color decomposition brings in additional powers of h which are not 
connected to the gradient expansion, but are related to the question to what extent 
color can be treated classically |15|. One should carefully differentiate between the 



limit of classical color and the classical limit in the kinetic sense which is defined 
through the gradient expansion. 

The transition back from QCD to QED is made by formally setting the color 
matrices equal to 1. Nonetheless the fields A^{x) are still operators in Fock space 
and neither commute with the Wigner operator nor with each other. The equations 



of motion for the Wigner operator are thus formally identical in QED and QCD, 
only that in QCD the (anti-) commutators refer additionally to the color indices. 



E. External field limit and Wigner function 

Equations of motion for the Wigner function (rather than the Wigner opera- 
tor) are obtained by taking ensemble expectation values of the operator equations. 
Writing for operators O working on W 

{[o, w]^) = [(6), {w)]^ + {[6, wu) - [(6), {w)]^ 

^[{d),{W)]± + C{d,W), (3.62) 

we can split the resulting equations into two parts: One part contains only one 
particle Wigner functions (W) and mean (gluon) fields (Vlasov term), while the 
other part describes the correlations C{0,X) (collision term). (For a discussion of 



the latter see the formulation with Green functions in [0.) In this subsection we will 
study the structure of the equations for the Wigner function under the assumption 
that the correlations can be neglected and that the fields F^^^ can be approximated 
by external fields {Ffj.u) = F^u- Then 

(qF^u) = {q)F^,, (3.63) 

and similarly for higher order operator products. This mean field approximation will 
then lead to a generalized Vlasov equation where collision terms, i.e. multi particle 
correlations, are neglected. 

For QCD the equations of motion for the Wigner function in an external gauge 
field are formally identical to those for the Wigner operator interacting with arbitrary 



gauge field operators, due to the non- commuting color structure. For QED the 
equations simplify somewhat since now the commutators of the gauge field with 
the Wigner function vanish, and only the non-trivial commutators with the Dirac 
matrices survive. 

In the external field limit the equations of motion for the electron Wigner function 
in QED read 

2m[l, W{x,p)]± = 2V^[Y, W{x,p)]i + tV^[Y, W{x,p)]^, (3.64) 
with (c.f. |8[) 

/■1/2 

'Pii = Pfj. + e ds sF^y{x + isdp) idp, (3.65) 

1/2 
rl/2 

T^n = dfj_ + ie J dsF^y{x + isdp) idp, (3.66) 

where e = —g is the electric coupling constant. Using Table ||, the two leading 
terms of the gradient expansion yield the semiclassical equations 

[ml -p^7'^,W(°)]± = (3.67) 

in zeroth order and 

[ml - p,7^ W^'^U + '£d,[Y, W^'^]^ = '£F,^ix)d;[r, W^'^]:, (3.68) 

in first order of the expansion. 

For QCD the zeroth order equation is again ( p.67|) , while the first order equation 
becomes after color decomposition 



singlet: [ml - Yp,, W^i'^]± + ^[7^ d,W^'^] 



'j^F^MY^W^'^^h, (3.69) 



octet: [ml - 7>^, W^^'^'']^ + ^[7^ d^W^''^'']^ 

+ ^r"'F',,d;[r, w^'^^u - '^r'^Kb', w^'^'h. (3.70) 

The equations for the color singlet distributions become identical to those of QED 
if we replace VT]") ^ VT^") and W^""^" N^{N^ - 1) VT^") with N^{N^ - 1) = 6 for 
QCD. 



The above mean field equations will be used in Sec. IV. 



F. Spinor decomposition 

We now return to the full operator equations and perform a spinor decomposition, 
using the following basis for the Clifford algebra: 

^, = l,^7^7^7^^^'^'^■ (3.71) 

In this basis the Wigner operator W is expanded as 

W{x,p) = ^ + t^'V + + Yl'A + ^^^'4-, (3.72) 

with the hermitian spinor components JF (scalar), V (pseudo-scalar), (vector), 
(axial vector) and S^y (tensor). These components are projected out from W by 



taking traces with the corresponding basis elements Fj [EH] . Inserting the expansion 



( p.72|) into Eqs. (|3.33|) and projecting onto the various spinor channels as described 



in ||23| one finds (after separating hermitean and antihermitean parts) 



2mJ^={n^,V^} (3.73a) 
2mP=[A^,>] (3.73b) 



and 



2mV^ = {n^,J^} + [A^4 



2mi^ = -[A^,P] + {^^44 



(3.73c) 
(3.73d) 
(3.73e) 



= [A„V^] 

o = {n^,>} 

= {n^, V4 - {n,, - £„^^.[a°, a"] 
o = [A^,j-]-{n^4j 



o = {^^,p} + [A^i: 



Here we defined the adjoint tensor tnv to Snv by 



(3.74a) 
(3.74b) 
(3.74c) 
(3.74d) 
(3.74e) 



(3.75) 



Up to the (anti-)commutator structure, these equations are identical to those derived 
in Ref. for QED in the external field limit. Note that by contraction with the 



Levi-Civita tensor Eqs. (|3.73| iii, |3.74| iii) can also be written as 



[A",v'']-{n,,ij + {n^,A} 



Eqs. ( |3.73|J3.7^ ) can be rewritten in the form of generalized transport equations 



and generalized mass shell conditions. To this end we introduce [30| a generalization 



of the Lorentz-covariant total time derivative (also known as the "drift operator" 
m dX{x,p)/dT = Pfj_d'^X{x,p) via 



{n^,[A^X]} or [A^{n^,X}] (3.76) 

and the generalized mass shell operator 

M^X = Arn^X - {n^, {n^ X}} + [A^, [A^ X]] (3.77) 

which will be used to generalize the classical mass shell conditions [m^ —p'^)X{x,p) = 
0. In the semiclassical limit and for Abelian gauge theories the generalized drift 
operator becomes the well-known Vlasov-operator p^d^ + ep'^dpF^^lx) while the 
generalized mass shell operator simply turns into — p"^. 

The details of the derivation of the generalized transport and mass shell equations 
are given in Appendix 0. One finds that only 8 of the 16 operators JF, V, V^, A^, S^^ 
are independent and satisfy generalized transport and mass shell equations while the 
other 8 functions are obtained from the solutions via simple constraints. For the set 
of independent functions one can choose either JF, "P, and S^y or and A^. We 
will give here the equations corresponding to the first choice; those for the second 



choice are found in Appendix D 3 . 



The complete set of covariant equations of motion for the spinor components of 
W{x,p) is: 

(a) Transport equations: 

[A^,{n^^}] = [A,A^,5n (3.78a) 

[A^,{n^P}] = KA,„i:n (3.78b) 

[A„, {n",4j]=[[p^ + n^,p, + n,],^] + e.^^^A^A^p] 

- ({4,, [A",p, + n,]} + {[p^ + n^, A°],4.}) (3.78c) 

(b) Mass shell equations: 



M^T={[p^ + A'^], 4J (3.79a) 

M2:P={[p^ + n^A'^],4j (3.79b) 

-[4,, [A", A,]] -[[A^, A"], 4.] 

-£a/3M.{[p" + n",A'3],P} (3.79c) 



(c) Constraints: 



0=[A^,J^]-{^^4J (3.80a) 
= {n^,P} + [A^4,]. (3.80b) 

(d) Defining equations for tlie dependent spinor components: 

2m4 = -[A^,P] + {^^4,} (3.81a) 
2mV^ = {n^, J^} + K>4,]. (3.81b) 

Tfie general structure tliat transport and mass sliell equations must be solved 
for only 8 of the 16 spinor components (either for JF, V, S^y or for V^, A^) while the 
other 8 components are obtained from constraints has apparently not been noticed 



before. It was not made manifest in the equations given in Ref. where transport 
and mass shell conditions were derived for the full Wigner operator but were not 
given explicitly for each spinor component. Compared to previous work Eqs. (|3.78 - 



3.81|) thus provide a considerable structural simplification of the covariant transport 



theory for spinor fields. 



G. Chiral Limit 



In this subsection we study the covariant kinetic equations in the chiral hmit of 
vanishing fermion mass. The resuhing equations of motion for the Wigner operator 
and its spinor components can be easily obtained by setting m = in Eq. ( p. 33 ) 



and in Eqs. ( |3.73| , p .7^ , ^.811 ). In particular one sees that now in Eqs. ( p.73| , p.73D 



Af^,V^ completely decouple from J^,V, and S^u- A closer look at the resulting 
"transport" and "mass shell" equations reveals, however, that they are no longer 
linearly independent but become linear combinations of each other. For m = it no 
longer makes sense to distinguish between transport and mass shell equations (which 
both contain the non-linear operators 11 and A in second order), and the transport 
theory is better formulated on the level of Eqs. ( p.73| , pT4D which are linear in the 
non-local operators. 

These equations can be further diagonalized by introducing a "helicity basis" via 



with \eQijkS^^^^^ = si^\ We find 



Oi 



= {n^,a'^W}, 
= [A^,a^W], 



(3.82a) 
(3.82b) 
(3.82c) 



(3.83a) 
(3.83b) 
(3.83c) 



as well as 



= [A^,/W]-{^^5W}, (3.84a) 
= {^„/W} + [A^sW]. (3.84b) 

This is the full set of covariant transport equations for gauge theories in the chiral 
limit, a^^-*, /'•^-' and s^^^ generate distribution functions for excitations with defined 
helicity/chirality. 

In a chirally symmetric theory both vector and axial vector currents should be 
conserved. Defining them by 



a/.(a;) = / j^Ai^,p), (3.85) 



one finds easily[| from Eq. ( p.83b|) the conservation laws 



= [D,ix),v^{x)], 

= [D^(x),a^(x)]. (3.86) 

The color singlet components of both currents (for which D^{x) reduces to df^) thus 
satisfy the usual conservation laws while the color octet vector and axial vector 
currents are covariantly conserved. 

In the recent literature chirally invariant transport equations were also de- 
rived in the framework of the Nambu-Jona-Lasinio model, in the external field limit 
and to leading order of the gradient expansion. The Lagrangian used there, 

>Cnjl = ti^Yd^^ + Glii^^y + , (3.87) 



^^Expanding in a Taylor series, all terms containing dp vanish after partial integration 
since V^^ vanish at infinite momentum. 



does not contain a coupling to gauge fields, but a self-coupling of the spinors which, 
if sufficiently strong, leads to spontaneous breaking of the chiral symmetry in the 
massless case. With the help of the auxiliary fields d = —2Gipil) and tt = —2Gipi'~^^il) 
the Lagrangian is rewritten as 



This Lagrangian is invariant under the chiral transformation 

^^^' = exp(-z75|)V' (3.89) 
0" ^ (t' = a cos X — sin X (3.90) 
n ^ n' = it cos X + cr sin X, (3.91) 

resulting in a conserved axial vector current. 



/^From this Lagrangian the authors of derived chirally symmetric transport 
equations in a vacuum state with dynamically broken chiral symmetry ((o"^ + tt^) ^ 
0). In contrast to our transport equations, these equations do not decouple 
and Afj, from JF, V and The coupling between the two sectors occurs via the 
vacuum expectation values a and vr of the auxiliary fields. However, according to 



Eq. ( |3.89| ), IT can always be made to vanish by an appropriate chiral rotation, with 



a compensationg chiral rotation on the spinor components of the Wigner density 
5T]. If one does this in the equations of Ref. |fT|, a takes over the same role as 



the bare mass m in our equations. In other words, the dynamical mass generated 
by spontaneous chiral symmetry breakdown couples the two sectors in exactly the 
same way as a non-zero bare mass m ^ 0. In the symmetric state, a = vr = 0, their 
equations decouple like ours. 



H. Classical limit 

Exact solutions for the kinetic equations in the free-field limit are given in Ap- 
pendix 1^. Here we study their classical limit, by expanding the spinor equations 
Eqs. (3.78-?^) to leading order in the gradient expansion. We will then consider 



explicitly the external field case for QED. 

With the generalized mass shell operator to zeroth order h, 

(M^) = 4(m2 - (3.92) 

which becomes just the classical mass shell operator, we find the following equations 



= 


(3.93a) 


iVfW'pW = 


(3.93b) 


M(°)'5(°j = 


(3.93c) 


= P^s^a 


(3.93d) 


= p'^p(o) 


(3.93e) 




(3.93f) 




(3.93g) 



The transport equations vanish identically, meaning that the equations to zeroth or- 
der do not contain dynamical information. From Eqs. (|3.93b|,P^93eD the pseudoscalar 



density proves to be zero. The equations are still not complete: they only give 
information on the functional p dependence of the spinor components. Only in the 
semi-classical approximation will we get transport equations for the spinor compo- 
nents and thus information on their x dependence. From the equations to first 



order h we find in addition mass shell conditions and constraints for the first order 
spinor components X^^\ that is information on their p dependence. 
To first order we find for the mass shell operator 

M2X = (M2)(°)xW + (M2)Wx(°) with 

= -'-lp^dp,[F^^{x),X^''>] (3.94) 



c 



and the generalized derivative 



first order 



2\p^D,{x),X^'^] - V9;{F.^, X(% (3.95) 



^Prom an explicit calculation we find further 



[[A^,A,],X] ^^'^^'■'^^'^O (3.96a) 



{b. + n„ A.], X} (3.96b) 
[W + n.,P. + n.], X] ^^^^*^^der ^^^^^^^^ ^(o)j (g^gg^^ 



This then leads to the following spinor equations in first order of the gradient ex- 
pansion 

\p,D^{x),P'^]-£p'^d;{F^,{x),P'^} = (3.97a) 



fjF,^{x),P''>] + I ({Fl{x),S^2} - {F,''{x),S^J}) (3.97c) 



(^2_p2)p(l) 



9 



9_ 

Ac 



(3.98a) 
(3.98b) 



(3.98c) 



SP'^^W = [/^,(x),^W] - |9'^{F.,(a:),^(°)} + |9,„{F"^(x), (3.99a) 
-2P,P(^) = [D^{x)t'^^] - |9,jF-^(x), 4°)} (3.99b) 



2mi(i) = 2p^£(i) - |9p.[F-'^(x), (3.100a) 

+ - |5,„{F"'^(x),5W}. (3.100b) 

Mass shell equations and constraints couple to the spinor components in first order, 
whereas the transport equations are equations for the zeroth order spinor compo- 
nents. Because of the vanishing pseudoscalar density V'^^^ the second transport equa- 
tion Eq. ( p.97| ii) becomes a constraint for u-^. Contracting Eq. ( |3.99| i) with and 



Eq. ( 3.99 ii) with p^e^"'^'^ and comparison with the transport equations Eqs. ( |3.97 i,iii 



we find as additional constraints 

[F^\Sf}]=Q (3.101a) 
[F'^'^,^^] =0 (3.101b) 

telling us especially that the field strength tensor commutes with the scalar density 

^(0)_ 



Thus we can conclude that the leading order in the gradient expansion is given 
through the mass shell and constraint equations Eqs. ( ^.9tj| ) and the transport equa- 
tions Eq. (3.97). (The additional constraints Eqs. ( 3.1Ul| i,ii) play no role to this 



order of the approximation.) To leading order the quarks (and electrons) can be 
interpreted as on-shell particles obeying a transport equation with gauge covariant 
drift and Vlasov terms, [p^D^{x), X] and g / {2c)pi^idpu{F'^'^ (x) , X} , respectively. The 
pseudoscalar density will contribute only to the next higher order. 

1. Classical solutions for QED with external fields 

For QED with external fields we get for the zeroth order formally the same 
equations as Eqs. ( |3.93| ), but now the spinor components are no longer operators. 
With -P^o) = this reads 

= (m^ - p2)jr(o) (3.102a) 

Q = {m'-p')Sf} (3.102b) 

= p^5W (3.102c) 

2mAf = ^p^Cf^ (3.102d) 

2mVf'^ = 2p^J^(°). (3.102e) 

To derive the first order equations we use that the commutators of the gauge fields 
and field strength tensor with the spinor components vanish, whereas the anti- 
commutators give a factor 2: 

0={p^d^ + ^F^uP^d;^ (3.103a) 
0=(p^d- + ^-F^pp'^d^) Sj^J + ^ (f/(x)5(°J + F/(x)5(°j) (3.103b) 



(m2-p2)^{i 



4c 



9_ 

Ac 



(3.103c) 
(3.103d) 
(3.103e) 
(3.103f) 
(3.103g) 
(3.103h) 
(3.103i) 



iFiom. the systematics of the expansion one finds that the hierarchy of equations can 
be solved successively for each order of the spinor components. The equations above 
for instance allow for an unambiguous solution for the classical spinor components 
X'^^\ without having to take into account higher spinor components. Therefore the 
hierarchy of equations can be truncated at any order for the spinor components. 
This remains even true for the more general case of QCD with external fields. The 
reason for this decoupling is simply that the operator to zeroth order % gives no 
contribution and therefore the equations to rfi^ order contain the rt"^ order spinor 
components only together with momenta p^. Allowing for correlations this remains 
true only if we replace multi-particle correlations by products of one-particle Wigner 
functions and external fields (Boltzmann Ansatz). 

To write down the equations for T^^^ , S^^J in their usual form with electric and 
magnetic fields we first use the separation into on-shell particle (5(po — E)) and 
anti-particle {6{po + E)) contributions for JF^^) and Sj^J according to Eqs. ( p.l02| i,ii) 
using the following identity 



5{p^ — rri' 



2bo| 



(5(po - ^) + 5(po + ^)) , ^ = Vp' + (3.104) 



In the following we will restrict ourselves to particles; the equations for anti-particles 
are given through po ~Po- The totally anti-symmetric field strength tensor is 



given in terms of electric and magnetic fields by 



/ 



—Ex —Ey —Ez 
Ex —Bz By 
Ey Bz —Bx 
E, —B„ Bx 



(3.105) 



Similarly we define the totally anti-symmetric tensor S 



I 



Sy (Jz 



Sx —(Jz CXy 

-ax 


with spin vector S and helicity cr. 

The dual tensor C'^'^ for S^'^ is then given by C^" 







-Sx Sy Sz 



\ 



Sz (^y (^x 



(3.106) 



i_^a(3nu^^^ in complete 



analogy to the dual field strength tensor J-'^^'^ = ^e^'^^^FaB P3]- ^^^^ is obtained 



from S^^ by exchanging S ^ cr and cr — *• — S. 

lYiom p^S^'^^'^^ = we find the constraints 

5^ = -vx o-W, V 
p ■ = 0. 



P. 

E' 



(3.107) 
(3.108) 



The first equation connects S^^^ and (t^^\ whereas the second equation is automat- 



ically fulfilled. 



Using 3-vectors we arrive after some algebra at the following equations for the 
dynamics of the system as governed by the scalar density and helicity cr'^^^ 
(coupling constant g = —e) 



[do + v + ^{E + v X B) ■ Vp]J^^°\x, p,E) = (3.109a) 
E[do + v-V + -iE + vx B)- Vp]crW(x,p, E) = 



--[E x{vx a^°\x,p,E)) + B X a^°\x,p,E)]. (3.109b) 

The other spinor components are dependent variables and connected to scalar and 
spin density through 

5(°)(x,p,E) = -vx cT^^\x,p,E) (3.110a) 

mA''Q\x,p) =p-(T^^\x,p,E) (3.110b) 

mA^°\x,p) = E(T^^\x,p,E) + px {v X a^°\x,p,E)) (3.110c) 

mV^ {x, p) = EJ^^^^ (x, p, E) (3.110d) 

mV^^\x,p) =pJ^^^\x,p,E). (3.110e) 

Eq. ( p.l09| ii) for the helicity or (with Eq. ( |3.110| i)) the spin density, is a generalization 
of the BMT equation . Eq. ( ^.109| i) for the scalar density corresponds to the 



Vlasov equation. The equations for the anti-particles are given through v — >■ —v. 
Eqs. (|3.109|j37ll0|) are still not sufficient to describe the theory, they have to be 



coupled to the classical Maxwell equations 



d.E'^" = f + j:^, (3.111) 
= 0. (3.112) 



jgxt is an external current and 



fix) ^ / d'p trYWix.p), f = jp^articics + J 



antiparticles ' 



(3.113) 



is the self-consistent current, coupling the field strength tensor in the Maxwell equa- 
tions to the Wigner function. We will get contributions only from the vector part 



and with the mass shell condition the current reads 



J 



u{0) ^ ± 



mJ (27r)3 



V 



(3.114) 



V (j^(o)(x,p,E) - J^W(x,p,E)^ 
Thus we find that 5^,^ couples only in next to leading order to the Maxwell equation. 
Therefore to zeroth order in the gradient expansion the self-consistent dynamics is 
governed by the scalar part JF(°) alone! 



2. Classical solutions for chiral QED with external fields 



Finally we study the massless mean field limit. The equations for .4.^(0) and 



V^(0) are now 



= - -F,,.d:] 
0= d^--F^,d; v^(°). 



(3.115a) 
(3.115b) 
(3.115c) 
(3.115d) 
(3.115e) 
(3.115f) 



Multiplying Eqs. ( p.ll5| iii,iv) with and using Eqs. ( |3.115| i,ii) we find that A^^^ 
and Vj^^ separate into a contribution from particles {5{po — \p\)) and anti-particles 



{^{po ~ IpI))- It further follows with Eqs. (|3.115| Lii,iv) that Af'^ and V^"'' can be 



,(0) 



expressed through A'^^ and Vq^^ 



A^'\x,p,±\p\) = e,A^^\x,p,±\p\) 
V^'\x,p,±\p\) = epVtf\x,p,±\p\). 



(3.116a) 
(3.116b) 



e„ denotes the unit vector in the direction of p. The transport equations in terms 



of the fields E, B, S,a are for particles (5(po ~ \p\)) 








a° + Cp ■ V + -{E + epX B)-Vp 



a° + ■ V + -{E + Cp X B) • 



A^o\x,P, \P\) 
vS'\x,p, \p\). 



(3.117) 
(3.118) 



For anti-particles one substitutes ip —ip and changes |p| — > —\p\ in the arguments 
of Aq'\Vq^\ The self-consistent current is now given by 



r(°)(x) 



d^p 2 
(27r)3p 



^ vl,'\x,p,\p\) + vS'\x,p,-\p\) ^ 



(3.119) 



yep {Vl,'\x,p,\p\)-v!,'\x,p,-\p\)) J 
Here A^j. will couple only in next-to-leading order to the self-consistent dynamics. 

The equations for V and S^j^u decouple completely from the equations for Afj, 
and V^, see Sec. [Ill G| . Therefore they will not couple to any order in % to the 
self-consistent dynamics so we do not have to treat them here. We should remark 
nevertheless that for JF, V and 5^,^ the hierarchy of the gradient expansion cannot 
be solved successively starting with the zeroth order components, but the equations 
of motion for the rfi^ component h couple to the n + 1^* component. 



IV. SINGLE-TIME FORMULATION 



We are now ready to make the transition to the single-time formulation of quan- 
tum transport theory, by taking energy moments of the covariant transport equations 
derived in the preceeding Section. As shown in Sec.|I|, a fully equivalent single-time 
description requires the calculation of all energy moments. 



A. Moment equations 



The starting point of our considerations will be again Eqs. (|3.27| , ^72^ ) for the 



Wigner operator. By multiplying from the left and right with 70 we can isolate the 
term {IIo, W} on the l.h.s. of the equations. Adding and subtracting then leads to 

2{no, W}=2m{j', W} - [7V, {n„ W}] - {7V, [^A„ W]} (4.1) 
0=2m[7°, W] - 2[zAo, W] - {7V, {H,, W}} - [7V, [zA„ W]]. (4.2) 

Through Ho Eq. ( [4.1D contains an extra factor po on the left, but no generalized time 



derivative [^Aq, W]. Eq. (|4.2| ) contains the generalized time derivative but no extra 
factor pq. Note that the generalized derivative A occurs only in commutators while 
n occurs only in anticommutators with the Wigner operator. 



For the spinor component equations (|3.73| , |3.7^) the analogous separation into 



equations with and without explicit po dependence yields 

{Ho, V°} = -2mJ^ + {Hi, V} (4.3a) 

{n°. A''} = 2m£°^- + ^"'^■'^[A,, V,] + {^^ (4.3b) 

{no,J^} = 2mVo-[A',5o.] (4.3c) 



and 



{n°, S^'^} = eo^Jk{2mAi + [A,, V]) - {^^ 5°^} + {n^ S''} (4.3d) 

{no,i°} = -{n„i^} (4.3e) 

{Ho, Vi} = {n,, Vo} + eoijk[A^, A"] (4.3f ) 

{n°, 5,0} = [A,, :F] - {W, 4} (4.3g) 

{Uo,P}^-[A\Coi] (4.3h) 



2mP = [A^, >] (4.4a) 

2m5o, = [Ao, Vi] - [A,, Vo] + £o^,■fc{^^ i'} (4.4b) 

2mVi = {ni,J^} + [A^5,,] (4.4c) 

2mA = -[Ao, :P] + {tf, jtoi} (4.4d) 

= [A^, VI (4.4e) 

= {n,, V,} - {n„ V,} - eoijk{[A', A'] - [A^ A']) (4.4f) 

= [Ao,J^]-{ff,5o,} (4.4g) 

= {n„ V} - [A°, £oi] + [A^ A,]- (4.4h) 

Before we can take po-moments of these equations we should first analyze the 
dependence of the non-local operators 11 and A on the partial derivatives dp^. We 
perform a Taylor expansion for the Schwinger string in dp^: 

MF.,(a; + tsdp) = ^ - a« + ^)| (qcd) (4.5) 



k=0 
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«=0 



For QED in the external field limit the higher order po derivatives thus couple 
simply to the higher order time derivatives of the external fields. As discussed 



in Sec. C|, for non-abelian gauge fields tlie above expansion generates a gauge 



covariant temporal gradient expansion of the transport equations. 

For the calculation of the energy moments the following identity is useful: 

{0 n < fc, 

(4.6) 

Here X stands for any spinor component of the covariant Wigner operator, and '"^X 
denotes the n**^ po-moment of X: 

Wx(a;,p)= r dpop^oXix,p). (4.7) 

J — oo 

Eq. ( [4.6|) is proved by partial integration, using that because of energy conservation 
we expect the Wigner operator to fall off at large po faster than any power. 

Whenever in the (anti-) commutators in the equations of motion one of the non- 
local operators 11 or A appears to the right of the Wigner operator, the momentum 
derivatives are defined as in ( ^.23| ) to act to the left in the sense of partial integration. 
For k momentum derivatives this introduces a factor (—1)'^. This motivates the 
definition of an alternating commutator through 

[A, B] [k even) , 
[A,B], = AB-{-1)''BA= I (4.8) 

[{A, 13} (A; odd). 
With these definitions one can work out the energy moments of the required 
commutators resp. anticommutators of the spinor components X with A and 11: 



J dpoP^[A^,Xix,p)] 



n 



NX(x,p)] + ^ [M(,)^(x,p),["-^']X(x,p)]^ (4.9) 



k=l 

n 



[Z}^(x),Wx(x,p)] + ^ [M(,)^(x,p),["-^']X(x,p)]^ , (4.10) 



A:=0 



dpoPo{'no,X{x,p)} 



k=0 



k+1 ' 



(4.11) 



dpop]^{Ui,X{x,p)} 



:{n,(x,p),NX(x,p)} + $: Q [%),(x,p),["-'=]X(x,p) 



k+l 



k=0 



k+l 



(4.12) 



(4.13) 



Here we defined the following non-local operators in 3-dimensional momentum space: 



^fiix, p)^P^i + ^ / ^^^^^ ^""^Fi^i^ + is^iSp) , 

A^(x,p)=;iL'^(x) + ^ /° ^ ^^s^H '"^^^^(^ + > 



(4.14) 
(4.15) 



M(fe)^(^,p)=^ f_^^^ds{-ihsf [aJJ) [^]F^,(a;o + yo,a; + i^sap) 



2/0=0 p 



ik 



2/0=0 



(4.16) 



2/0=0 P 



2/0=0 



(4.17) 



Note that we have used for the 3-dimensional momentum space versions of the 
operators 11 and A the same symbols as for the 4-dimensional ones; which one is 
meant in a particular equation should be obvious from the context. 

For a given power n of the term Ho gives rise to moments which are one order 
higher than those from all other terms. This term thus indeed couples different 
orders of the moment hierarchy so that it will not trivially truncate. 

The moment equations are now given by 



k=l 



+{7V,[^M(,)„["-'=]#],}) 
(constraints) 



(4.18) 



+ [7V,[^M(,)„['^-'=W],]). 
(dynamics) (4.19) 

The first of tlie two sets of equations, which couples the n lowest moments to the 
n + 1*^* one, contains only derivatives with respect to the spatial components of x 
and no time derivatives. We can formally write it as 

n 

["+ilH^(a.,p,xo) = E Ol^{x,p,x,- d,, d,) m{x,p,x,) , (4.20) 

fc=0 

which shows that it forms a set of constraints. 

The second set of equations determines the time derivative of the n**^ moment in 
terms of all lower moments including the rfi^ one, which we write formally as 

d, M#(a^, p, xo)+Q^,)(x, p, xo; 9., d,) M#(a^, p, x,) 

n-l 

=- E Q(fc)(x,p,xo;a,,a,) ^'W{x,p,xo). (4.21) 

k=0 

This is a time-dependent partial differential equation for the ra*^ moment, with a 
source term from the lower moments, which we interpret as a single-time transport 
equation. The equation for n = is special because in this case the source term 
on the r.h.s. of (|4.21| ) vanishes. The operators C^(fc), Q^k) functionals of the 



gauge fields and contain arbitrary powers of the spatial coordinate and momentum 
derivatives. 



The spinor decomposition of the moment equations is given in Appendix ^ 



B. BGR equations 

Limiting ourselves to external fields and considering only the lowest (zeroth) 
moments we can reduce the dynamical equations ( F9| - |F16| ) for the spinor components 



to the BGR equations p4| , |38| . For QED with external fields we find for the non-local 
operators 

/external n / \ 

dpop-, [A,,X{x,p)] ^ V, Nx(x,p) + Y: [f^)K(k), (4.22) 

/ dpoPn^o,Xix,p)} m 2["+i]X(x,p) + 2^ pL(,)o (4.23) 

J dpopnU,,X{x,p)rm2V. Wx(x,p) + 2f: (^^^L(fc), ["-'^]X(x,p) (4.24) 

with the definitions p.65| , ^.6^ ) for V^, and the following operators and 

1/2 

Ki^k)^.=g I ds{-isf[df^F^,{x + isdp, xo)di + '-d^^-^^ F^^i^x + isd^, xo)] , (4.25) 

-1/2 ^ 
1/2 

Lik)^.=g [ ds{-isf^\df;jF,^{x + isdp, xo)dl + -5^-1)^0^(2; + isd^, xq)]. (4.26) 

-1/2 

The BGR equations are thus given as follows 

= 2m[7°, - 2iPo - 2{7°y , Pi '"IW^} - [7°^ , zP^ (4.27) 

In spinor decomposition they are 

2m Mp=P^ [°U^ (4.28a) 

2m M5o.=I)o - [°lVo + 2e,i,kV' [oU'^ (4.28b) 



2m ^^^Vi=2V^ [°]J^ + V (4.28c) 

2m [°Uo=-Po + 2V' [°l£o* (4.28d) 

0=p^ My^ (4.28e) 

0=2P, - 2Vj !°lV, - £oiifc(^?° - V [°U°) (4.28f) 

0=1)0 - 2P^ [°l5oi (4.28g) 

Q^2V^ - P'' [°l£oi + 'D^ (4.28h) 

The BGR equations therefore are just the dynamical equations for the zeroth mo- 
ments for the case of QED with external fields. The connection with the spinor 



components fi,gi, i = 0, 1, 2, 3, as introduced by BGR in [0] is 



/o = '°lVo, /i = -[°Uo, /2 = [°lP, /3 = ™^, (4.29) 



go = -[°U, gi = -™V, g2 = ^'^S, gs = -^'^cT. (4.30) 



The hierarchy of moment equations Eqs. ( [4.18| , ^?T9D will in general not truncate. 



The initial problem, namely that the covariant Wigner operator involves the fermion 
fields at all times and thus cannot be properly initialized at t = — oo, has now been 
reformulated into the task of solving simultaneously infinitely many equations as 
initial value problems. On the other hand the BGR equations (|4.27| , |4?28| ) can not 



describe the theory completely since they lack information on all higher moments. 
How we can truncate the hierarchy naturally and what role the BGR equations 
play with respect to the infinite hierarchy of moment equations will become clear in 
the next subsection where we discuss the moment hierarchy in the limit of external 
fields. 



C. Moment hierarchy for external fields 



Since the Dirac equation and its adjoint for (time dependent) external fields are 
inhomogeneous first order differential equations for the wave functions, we need to 
specify as initial conditions only the fields ip, ip'^ for all x at some initial time xq = ti. 
In terms of moments it should be sufficient to specify '^'Vr at some initial time t,. 
From the structure of Eqs. ( ^4.20|) it is obvious that there exists a recursion relation 
which allows to express all higher order moments ["lW^(cc,p, Xq), n > 1, through the 
lowest moment MW^(a;,p, xq). Let us write it schematically as 



^W{x,p,xo) = Vn{x,p,XQ]d^,dp)^^^W{x,p,XQ) , n > 1, 



(4.31) 



where Vn is a combination of products of the operators O";^.-) and thus contains 
derivatives with respect to x and p but no time derivatives. With the help of ( |4.31|) 
one can rewrite (|4.21|) for n > 1 as 



do {Vn ™ W^) = - E Qh'^k ^'^W, ,n>l, 



(4.32) 



A;=0 



or 



k=0 



W, ,n>l. 



(4.33) 



Ehminating on the l.h.s. of (|OB| ) do with the help of Eq. (|]^) forn = we 
obtain 



fc=0 



W, ,n>l. 



(4.34) 



This is, for each n > 1, a new constraint equation for ^^^W. 

We have thus replaced the set ( [4.21|) of single-time transport equations by a 
single transport equation for (the case n = in ( [4.21|) ) plus an infinite set of 



new constraint equations ( [4.34| ) for n > 1. The single remaining transport equation 
reads explicitly 

2[zAo, = 2m[7°, 1°%] - {7V, {n„ - [7V, [^A^, . (4.35) 

This is the non-abelian generalization of the BGR equation ( [4.27|) for QED. 



Eq. (|4.35|) is a first order differential equation in xq for ^^W which can be solved 
uniquely by specifying ^^^W{x,p, xq) at xq = ti. But what is the meaning of the 
additional constraints ( [4.34|) ? 



We will now show that they are redundant, i.e. turn into identities after using 
the constraints ( [4.31| ). A brute force proof of this statement was given in for 
the case of a simple scalar field theory. For QED the statement was proven only 
for the lowest non-trivial values of n. The reason for this was that for higher n the 
explicit calculation of the operators P„ in ( [4.31| ) becomes extremely cumbersome. 



We will here give a more elegant general proof which is also directly applicable to 
the non-Abelian case which exploits the "equivalence" (with respect to taking higher 
Po-nionients) of the two covariant equations from which the single-time transport 
and constraint equations were derived. 

The proof is made more presentable by introducing the following notation for 
multiple commutators: 

A_X=[A,X], A+X = {A,X}, A.B.X =[A,[B,X]] etc., (4.36) 

where the commutators or anti-commutators extend to everything to the right of 
the operator A±. For the spinor structure we will use in addition 

g' ^ 7°, ^ 7V (4.37) 



Let us now look at the integrands under the moment integral / dpo Po ^he two 
equations which gave rise to Eqs. ( [4.20| , ^.2lD : 







W, 

w 



(4.38) 
(4.39) 



We will now consider in addition Eq. ( |4.39|) multiplied with 2po (which can be written 
aspo+[Eq. (IH] ) 



2m^° - 2iAo_ - g\Iii+ - i glA^ 



(4.40) 



After taking moments / dpoPQ, Eqs. ( [4.39| , ^l40D give dynamical equations for the 
moments n, n + 1 while Eq. ([4. 38]) gives constraints which tell us how to eliminate 



,n. 



the moment n + 1 through the moments k = 0, 

The idea of the proof is to push in Eq. ( [4.4(JD the operator po+ to the right until 
it stands in front of W, then to use Eq. (|4.38| ) to ehminate po+W and afterwards 
exploit Eq. ( [4.39|) to eliminate the dynamical part Aq_W. The result will be an 
identity. On the level of moments this is equivalent to the elimination procedure 
discussed at the beginning of this subsection. 



We first push po+ in Eq. ( [4.40| ) to the right and insert Eq. ( ^4.38|) . Noting that 



g!^ commutes with po+, ^/x-, we find 



mg_ 



lAo- - lglU,+ - y_A,_{2po+W) 



+ 



2(iAo-Po+ - Po+'^Ao-) + gX(Jli+po+ - po+n^H 



w 



mg''^ - 2Ao_ - hlUi+ - \glA,_ 2mgl - 2N^+ - glY{i+ - g'A,^ 



W 



W. 



(4.41) 



Now (i.e. after having inserted ( |4.38| )!) we push Aq- again back to the right and 
use Eq. for Aq.W: 



2mgl - 2No+ - giUi+ - g^A^ 



mg 



29+^^1+ - 29- 



W 
W 



+^V(ni+Po+ - Po+n»+ + iAo_iAi_ - iAi_iAo_^W 
+gi(iAi^Po+ -po+iAi^ + iAo-U,+ - Ui+iAo^W. 



Expanding and combining terms with the same spinor structure leads to 



2m\g^^gl - glg"_) - 2mN,^{g'_ - g' ' 
-mU,^{g''_gl + g\gl - glg^_ - glg\) 
-miA,_{g^_g\ + g^_gl - glgl - g\g'_) 
+^;[ni+no+ - + ^Ao-iA,_ - iA,_iAo- 

+c/i[zAi_no+ - no+iAi_ + 'iAo-ni+ - Hi+iAo- 



+W-9i - g\gl)[iA,-iA,. - n,+n,+] 



-i+J 



W. 



(4.42) 



(4.43) 



As can be easily checked the following relations for the spinor matrices hold: 



i9'L9l - 9WW = 



(g'.gl - gW. + glgl - gXg'.)W = 



(4.44) 
(4.45) 



{g-g- - glg\ + g\gl - g-g-)w = o (4.46) 

{g\gi - glg'-)W = 2{g^Wg^ + gmg') (4.47) 
{g'-gi - gy-)W = 2{g'Wg^ - g^Wg'). (4.48) 

Therefore the terms containing the mass m will vanish. The rest can be combined, 
using the explicit form of the 7 matrices, to give 

= [(n^+ + «A^_)(n,+ - «A,_) - (n,+ - «A,_)(n^+ + zA^„)] 7Vw^7V- (4.49) 

With the help of Eqs. ( P7|) we finally arrive at 

-7°7'^iy7V[p. + - ^K,P^. + + ^A^]. (4.50) 

In Appendix ^ we show that the commutator vanishes: 

[p^ + U^ + iA^,p, + U,-iA,] = 0. (4.51) 

Eq. ([4.50|) is therefore an identity. 

In the case of external fields we thus have a natural truncation of the transport 
hierarchy on the level of ^'^^W. Its time evolution is uniquely determined by Eq. ( [4.35| ) 



once W{x,p,t)\t=to has been specified. All higher moments can be found through 
successive insertion of the solution into Eq. ( ^.181 ). These "constraints" are differen- 



tial equations in x and p without time evolution. For the external field problem the 
"initial value problem" for the quark Wigner operator has thus been solved, and we 
have also given an explicit prescription how to find the higher moments. The latter 
was missing in the work of Refs. p4|-p^, although the BGR equations Eqs. ( [4.27|) 



properly describe the dynamics of the zeroth moment. Higher moments do, however. 



contain important physical information even for external fields. For example, the 
energy momentum tensor is given by the first moments (for explicit expressions see 
||27| ]) and can therefore be calculated only with the help of the constraint Eq. (4.18) 
once the solution for the zeroth moment has been obtained. 



D. QED and external fields 

In this subsection we make the results obtained above explicit for the case of 
an ahelian external gauge field, i.e. for QED with external fields. We will show 
the explicit connection between the covariant theory and the moment equations in 
leading order of the gradient expansion and demonstrate that the correct Vlasov 
(including spin precession terms) limit is recovered from the single-time transport 
equations. 

To zeroth order in the gradient expansion the equations for the moments are 



[n+l]yO(0) 


= m 




(4.52a) 


["+i]jf(0) 


= m 


Nyo(o) 


(4.52b) 


[n+l]yj(0) 


= P' 


NyO(o) 


(4.52c) 




= P' 


N jf(o) 


(4.52d) 


pi Myi(o) 


= 


[n]yi(0) 


(4.52e) 


[n+l]p(0) 


= 




(4.52f) 


Hp(o) 


= 




(4.52g) 


pi W-p(O) 


= 




(4.52h) 



and 



[n+ll^fc(0) ^ _^ [n]£Ofc(0) j^^k [n]^0(0) 

l]gjm ^ ^^Oijfc N^(0) _pk [n]^Oi(0) [n]^Ofc(0) 



n+ 



[n+l]^0(0) ^ N^i{0) 

[n+l]^Oi{0) ^ [n]^ii(0) 

^ N^O(O) ^ ^. [n]£Oi(0) 

0=p, W50i(o)_ 



(4.53a) 
(4.53b) 
(4.53c) 
(4.53d) 
(4.53e) 
(4.53f) 
(4.53g) 



Eqs. ( [4.52| vi,vii,viii) show that vanishes for all n identically. In addition 

we get that V^°-* can be expressed by JF^^) and ^[f-* by S^^} . Eqs. ( |4.52| ) can be 
solved now for MjF(o) and MV^^) by specifying the lowest moments [°1jf(°) and Mv°(°) 
{E = y/p'^ + m'^ > 0, m > 0): 



[2n]y0(0) ^ ^2n [0]yO(0) [2n+l]y0(0) ^ ^^n+l ^ [0] jp(0) 

m 



[2n]y(0) ^ ^2nP_ [0] jr{0) [2n+l]y(0) ^ ^2n+lP [0]yO{0) 

E 



m 



(4.54) 



[2n] -p(O) ^ _g2n [0] j-(0) 



Using 



[2n+l]jr(0) ^ _g2n+l^ [0]yO(0) 

E 



^ [0]yO(0) ^ [1] -p{0) 



(4.55) 



we can also write 



[2n]y0(0) ^ _g2nj_ [1] jr{0) [2n+l]y0{0) ^ _g2n+l:^ [0] j^(0) 



[2n]y(0) ^ ^2nP_ [0] jp(0) [2n+l]y(0) ^ ^2n+l^ [1] jp(0) 



(4.56) 



[2n]j-(0) ^ ^2n [0] jp(0) [2n+l]j-(0) ^ ^2n+l^ [1]jf(0)^ 

thus specifying rather MjF(°) and 1^1 JF^'') instead of MjF(°) and My^'-^''. The connection 
with particle (upper sign) and anti-particle (lower sign) distribution functions from 
the Lorentz-covariant theory is found to be 

^(x,p, ±E) = E [01^(0) ± m [°lV°(°) = E ± (4.57) 

Thus the particle and anti-particle contributions of the covariant theory J-'{x, p, ±-E) 
translate into the sum and difference respectively of ''^^JF and '^^JF in single-time 
formulation. This is immediately obvious from the fact that is derived from 
as an integral over all energies, that is particle and anti-particle contributions 



'ljF(x,p)= J dpoJ^{x,p), (4.5^ 



whereas '^'jF has an additional weight po, leading to the difference between particles 
with positive and negative energies. 

The transport equations for the zeroth order moments in the semi-classical ap- 
proximation are now given by (see Appendix ^ 



0=E'' 



E{d' + -V„ ■ ^;)W + (p . V + -(p X S) ■ Vp)(^[°i^(°)) 

c c , 

E(9° + -Vp ■ ^;)(eM^(°)) + (p . V + -(p X s) ■ Vp)W^(°) 



(4.59) 
(4.60) 



Adding and subtracting gives us back exactly the transport equations we get from 
putting the distribution function of the covariant theory on the mass shell 



inserting into the transport equations of the covariant theory Eqs. ( p.l03| i, |3.109| i) and 



integrating over jdpo p^. This shows by example the equivalence between covariant 
and single-time formulation. 

Similar we get for W^Co) and '"^iS^^-* the following equations, if we specify the 
lowest moments of W5^^(°), that is [°1(t^(°) and 

[2n]^0(0) ^ ^2nlp . [0]^(0) ^Qla) 

m 

[2n]^(0) ^ ^2nl [1]^(0) + p x (p X (4.61b) 

m \ J 



[2n] giO) 


= E^^p X - 


[2n]^(0) 


^ ^2n [0]^{0) 


[2n+l]^0(0) 


m 


[2n+l]^(0) 


^ J^2n+l}_ ( 

m \ 



(4.61c) 
(4.61d) 

WW 

(4.61e) 

_+px(px-^)j (4.61f) 

[2„+l]^(0) ^^2n+l_^_ (4.61h) 

As before we have for the connection with the particle and anti-particle distribution 
of the covariant theory 

o-(x,p, ±E) = E Mo-W ± (4.62) 

The transport equations as derived in Appendix |^ are 



= + -Vp ■ E)—^ + (p . V + -(p X B) • Vp)[°W(°) 

C c 



+ ^ p-(£;-^) + Sx[V(°)) (4.63a) 



p p [i]/t(o) 



(4.63b) 



Adding and subtracting gives us again the transport equations for the particle and 
anti-particle distributions of the covariant theory Eqs. ( p.l03| ii, |3.109| ii) after integrat- 



ing over / dpo Po and employing the mass shell condition cr = cr(x, p, ±£')(5(po T -£')• 
Eqs. ( |4.63| ) are therefore again the generalized BMT-equations Eqs. (|3.109| ), but now 
formulated with moments. 

In the case of QED with external fields we thus showed the explicit connection 
between the formulation with moments and the covariant theory. To this end we had 
to solve the moment hierarchy exactly in the classical limit. This solution proved 
then to correspond to the classical mass shell condition of the covariant theory. 



V. CONCLUSIONS 

In this paper we have established the connection between the covariant and the 
so-called single-time formulations of the quantum kinetic theory, which is based on 
the dynamics of the quantum mechanical Wigner operator. We discussed the role 
of the time parameter in the density matrix, and we could show that a formulation 
with a single time parameter necessarily leads to a hierarchy of equations for the 
energy moments for the Wigner operator. For these moments we found a direct 
connection with the time derivatives of the wave functions. 



A covariant field theory is equivalent to the complete set of its po-moments. 
The impossibility to formulate the covariant transport theory as an initial value 
problem, which we discussed in Sec. ||, thus translates into the impossible task of 
solving an infinite hierarchy of moment equations. The hierarchy structure lends 
itself, however, to a systematic approximation scheme via successive truncation. 
For fermion dynamics in external classical gauge fields the truncation problem can 
be solved exactly, in the case of Dirac fields at the level of the lowest po-HioHie^t, 
the single-time Wigner function W3{x,p,t). This reflects the fact that for external 
gauge fields the fermionic theory is a single-particle theory whose dynamics can be 
solved uniquely once initial conditions for the fermion wave functions iP{x,xq) at 
Xq = ti have been specified. 

To find the equations of motion for the moments we first reviewed the covariant 
transport theory, extending some of the results obtained previously in [ 
introduced a compact notation for generalized non-local momentum and derivative 
operators which allowed us to cast the equations of motion for the Wigner operator 
into a form which permitted further manipulation. We then performed a gradient 
expansion followed by a color decomposition and established the connection be- 
tween the (local) momentum operator and covariant derivative with their non-local 
generalizations. 

To find the proper mass shell and transport equations for the Wigner operator 
we had also to perform a spinor decomposition. Thereby we could eliminate half 
of the initial degrees of freedom to end up with a system of equations for only 
eight independent spinor components, whereas the other eight components could be 
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expressed as dependent variables. This improves on previous covariant treatments 
were the system of mass shell and transport equations always contained redundant 
information. We derived this system of mass shell and transport equations plus 
constraints on the one hand for scalar, pseudo-scalar and tensor components, which 
are a suitable basis of spinor components for massive quarks; alternatively we also 
derived mass shell, transport and constraint equations for vector and axial vector 
components, which are more suitable as independent components for the case of 
vanishing bare quark masses (chiral limit). 

That our equations have the proper classical limit was proven by an explicit 
gradient expansion. We could also find the expected vacuum solutions where the 
pseudo-scalar and spin density did vanish identically. For chiral theory we found 
the expected conservation of vector and axial vector current. We also compared our 
results with recent investigations of quantum transport equations for the Nambu- 
Jona-Lasinio theory. We showed that the failure of those equations to decouple 
scalar, pseudoscalar and tensor densities from the vector and axial vector ones was 
related to the spontaneous breaking of chiral symmetry in that theory, with the 
dynamically generated mass acting like a non-zero bare mass in our case, thereby 
coupling the two sectors. 

After the discussion of the covariant theory we considered the expansion of our 
equations in energy moments. This moment expansion was performed with the 
help of a partial gradient expansion in d^^ and resulted in an infinite hierarchy 
of dynamical equations for the moments, which are coupled by constraints. The 
equivalence of both formulations was given by construction. For external Abelian 



fields to leading order of the gradient expansion we showed the explicit connection 
of the particle and anti-particle distributions of covariant theory with the zeroth and 
first moments in the equal time formulation. Finally we were able to show that for 
the external field case we need only the dynamical equations for the zeroth moments. 
All higher moments can then be successively calculated from the lowest moments 
via the constraint equations. This leads for the external field case to a natural 
truncation condition for the moment hierarchy where the minimal truncation, that 
is ignoring all but the lowest moments, just corresponds to the BGR equations. 

The question now arises whether now the equations of motion for the single- 
time Wigner operators can be finally applied for the (numerical) description of a 
QGP or electron-ion plasma. For the transport equations with the full operator 
structure the answer is certainly no in practice. Further work is required to analyze 
correlations and collision terms in this case which should be obtained by taking 
an ensemble average and using appropriate factorization of expectation values. For 
external fields, however, where the equations reduce to transport equations without 
a collision term, quantum transport theory has now been formulated in the single- 
time approach in a way which allows for immediate numerical implementation, for 
both QED and QCD. For each color and spin component one transport equation 
must be solved, together (but not simultaneously!) with as many time-independent 
constraints as given by the highest order of po-moments which one wants to know. 

Quantum transport theory beyond the external field limit and the extraction of 
suitable collision terms remains a great challenge. For gauge theories the Wigner 
formalism, through the occurrence of Schwinger strings, strongly suggests to use the 



radial gauge |^ for which perturbation theory has only recently become viable after 



a set of consistent Feynman-rules was established in Ref. E^]. Infrared divergences 



and the need for resummation in non-abelian gauge theories with massless degrees 
of freedom have, however, so far hampered practical progress. 

APPENDIX A: COLOR MATRICES 

The generators of SU{3) in the fundamental representation, ta = |Aa (where 
are the Gell-Mann matrices), satisfy the following identities 

= e (Al) 

{e, t^} = — + hd^'^H" (A3) 
3 



1 



= ^ + /irt'^j . (A4) 

jabc ^abc ^-j^Q totally auti-symmctric and totally symmetric structure con- 
stants of SU (3). 

The traces of the generators are given by 



tr t'^ = (A5) 

tr tH^ = —6''^ (A6) 

tr eth" = ^ (d"''" + if^') = ^h''^" (AT) 

tr tHh^f = -^5^^- (A8) 



One also has 



r^^ = -i^tr (r,^^) (A9) 

d"'''' = p tr {{f, t^'ji") (AlO) 

f"^^ = (All) 

d"^^ = (A12) 



APPENDIX B: GAUGE TRANSFORMATIONS 

Let S{x) be a gauge transformation in U{1) or SU{3): 



(Bl) 



(For U{1) there is only one generator, to = 1.) Then the following transformation 
rules apply: 



ip{x) - 


-> S{x)ip{x) 




(B2) 


iP^ix) - 


-^7jj''{x)S-\x) 




(B3) 


M^) - 


-.Six) 




S-\x) 


(B4) 


D,{x) - 


S{x)D^{x)S-\x) 




(B5) 


U{x,y) - 


^ S{x)U{x,y)S-\y) 




(B6) 


g{x + y/2,x-y/2) - 


^S(x)Q(x + y/2,x-y/2)S-\x) 


(B7) 


W{x,p) - 


S{x)W{x,p)S-\x) 




(B8) 




S{xy^^F,,{z{s))S-\x) 




(B9) 


[D^ix),W{x,p)]- 


S{x)[D^ix),Wix,p)]S-\x) 




(BIO) 



^ S{x)U^S-\x) (Bll) 

{U^,W{x,p)} ^ S{x){U^,W{x,p)}S-\x) (B12) 

^ S{x)A^S-\x) (B13) 

[A^,#(x,p)] ^ 5(a;)[A^,#(x,p)]5-i(x). (B14) 



APPENDIX C: LINK OPERATORS 



We recall the formula for the variation of link operators with respect to their end 
points p^: 



6U{b, a) = 6b igA{b)U{b, a) - igU{b, a)A{a) 6a (CI) 
-tg ds U{b, z{s)) F^^{z{s)) U{z{s), a) {b - {6a + s{6b - 5a))^ 



z{s) = a + {b — a)s, 



This yields the following explicit expressions: 



djf-^U{x+,x)=—-igU{x^,x) — 2[/(x, x+)A^(x+)f/(x+, x) 



1/2 



+ / dsil + 2s)^^^F,^ix + sy)y'- 



d+yU{x,x-)-- 
d-^U{x+,x)-- 
d-^U{x,x-)-- 



--Ug [~A^{x) + f ds (1 + 2s) + sy)y^] U{x, x. 

Zi J — '1/2 

--\igU x) \a^{x) + [^'^ ds (1 - 2s)[^]F^^(x + sy)y'' 

1 



-A^{x)^2U{x,x^)A^{x^)U{x^,x) 



+ /° ds{l-2spF,^{x + sy)y'']u{x,x.). 

^1/2 J 



(C2) 
(C3) 
(C4) 



(C5) 



For the derivative of the Schwinger string ^^^Fi,^{x + sy) (see Eq. ( pj.24D ) with respect 
to y one finds 



dyJ^'^Ff.^ix + y) = U{x,x + y) [dya - igAa{x + y), F^i,{x + y)]U{x + y,x) 



-ig dssy^ + sy), [^'^^^(x + y) 



(C6) 



Its covariant derivative with respect to x is 



V^{xpF^p{x + y)=[D,{x), + y)] 



U{x,x + y) [V^,{x + y)Fai3{x + y)]U{x + y,x) 



-ig dsy-" y^\F,,{x + sy)^^^F^p{x + y) 
Jo L 



(C7) 



Combining the last two equations one finds 



V,ixpF^p{x + y) = 

dy^^^^F^pix + y)-ig f' ds (1 - s) y^ Y'F^^ix + sy), + y)\ . (C8) 



APPENDIX D: SPINOR EQUATIONS 

In this Appendix we combine the spinor equations ( p.73| , p.7^ into generalized 
mass shell and transport equations for the spinor components. Multiplying the 
first three of equations in ( |3.73| ) and ( |3.74| ) each by 2m and inserting the last two 
equations ( |3.73| ) yields 



4m2^ = {n^, {n^^}} + {n^, [a,,^^} 

= -[A^, [A^P]] + le^''''^[A„ {n,,5,4] 



(Dl) 
(D2) 



4771^4, = [A^, {n,,^}] - [A„ {n^,^}] + [A^, [A^Jua]] - [A,, [A",4„]] - 



(D3) 



as well as 



= [A^, {n^ T}] + [A„ [A,, 5^] (D4) 

= {n^, [A^p]} - ^e^^^^^in^, (D5) 
= {n^, {n,, j^}} - {n,, {n^, j^}} + {n^, [a",5,«]} - {n,, [a°,4„]} + 

£„^^4A-, [A^ V]] - ^e^p^^e^^'^lA'^, {n^, Ss.}]. (D6) 

With some further manipulations the first three equations will give rise to generalized 
mass shell constraints, the last three to generalized transport equations for JF, V, and 
Sfj^iy. In terms of the solutions to these equations the remaining spinor components 
Afj, and are given by the last two equations ( p.73| ). 



For the further manipulations the following identities will be useful: 

[A^, [A., X]] - [A., [A^, X]] = [[A^, A,], X] (D7a) 

{n^, [A., X]} - [A,, {n^, X}] = {[p, + n^, a.], x} (D7b) 
{n^, {n,, X}} - {n,, {n^, x}} = [[p, + u„ p, + u,],x]. (D7c) 

The additional factors ~ p^ on the r.h.s. arise from the fact that the commutators 
in these equations refer to the color and Fock space structure of the operators only 
and do not consistently include also the momentum derivative operators which, 
according to (|3.23| ), were defined to act only to the left when standing to the right of 



a Wigner function. Splitting the generalized momentum operator (p.25|) according to 



= Pn + TTfj, where the second term vr^ (like A^ in ( p.26|) ) contains only momentum 
derivatives, but no momentum factors, Eq. (P7b| ), for example, therefore has to be 
evaluated as follows: 

{p^ + TT^, [A„ X]} - [A„ {p^ + TT^, X}] 

=+p^(A,X) -p^(XA,) + {A,X)p^ - (XA,)p^ 



+{7r^,[A,,X]}-[A,,{7r^,X}] 



=2[p^, A,]X - 2X[A,,p^] + [vr^, A,]X - X[A 



+ n^, Ajx - x[A„p^ + n^] 

= [[p^ + n^,AJ,X], (D8) 



and similarly for Eq. (|D7c| ). In the intermediate expressions the round brackets 



indicate, where ambiguous, how the momentum derivatives act via ( |3.23|) ; since the 



momentum derivatives in vr^ and A^ commute the remaining terms are unambiguous. 
We will now use the generalized quadratic mass shell operator 

M^X = 4m'X - {n^, {n^ X}} + [A^, [A^ X]] (D9) 

to combine equal spinor components on the left hand sides of the equations. 

1. Mass shell equations 

We now derive the mass shell constraints ( |3.79| ). To get the mass shell condition 
for JF with the operator from above we bring the term {11^, {11^, JF}} in Eq. ( P9|) 
to the left and add on both sides [A^, [A^, JF]]. Furthermore we employ the constraint 
Eq. ( |3.74| iv) together with the anti-symmetry of iS^^/ in the form 

[A^, [A^^]] + {n,„ [A^^ni = [A;., {n,,^'^'^}] - {n„ [a^,^^} , (dio) 

and then use ( |D7b| ). 



The mass shell condition for V we obtain by subtracting {n^,{n^,P}} from 
Eq. (021) and using 



[A„ {n., tn] - {n„ {n^ p}} = [a^, {n,, - {n., [a^, &^]} , {mi] 



which follows from Eq. ( p.74| v). We then again use ( D7b ). 



For the mass shell condition for finally, we rewrite £Q,/3^,^{n", [A^jP]} as 
e„^^,{n°, [A^,V]} = e^p,,{[A(', {n",P}] + {H", [A^,V]} - [A^ {U",V}]), (D12) 
and use the constraint ( ^.74| v): 

e^^^^n", [A^,V]}=e^p,. ({n°, [A^,V]} - [A^ {U",V}]] 



-[A^ [A,, 4^]] - [A^ [A,J^,]] + [A^ [A^J,.]]. (D13) 

One further employs 

[A^, {n,, j^}] - [A,, {n^,^}] = [A,„ {n.,^}] - {n., [a^, j^-]} 

-([A,,{n^,^}]-{n,„[A,,^]}; 
+{n,, {n", 4j} + {n^, {n", 4.}} (dm) 



where the constraint ( 3.74| iv) was used in the last line. Inserting now Eqs. (pi3| , [nT4[ ) 



into (^) and using Eqs. (^) leads to the mass shell condition for 



2. Transport equations 

In the following we want to express the transport equations with the generaliza- 
tion 

[A^,{n^X}] (D15) 

of the Lorentz covariant (proper time) derivative fn-^X = p^d^X (r: proper time, 
Xq: global time). For JF Eq. (p4|) is already in the desired form. The transport 



equation for V follows by acting with on Eq. ( |3.73| v). For S^,^ we use the 
following identity 

{n^, [A",5„,]} + {n,, [A",4„]} = {n^, [A",5„,]} - [A", {n^,5„,]} 

{n,,[A°,4,]}- [A",{n,,4„]}. (D16) 



3. Equations for the vector and axial vector densities 

Alternatively one can derive transport equations, mass shell conditions, and con- 
straints for the eight vector and axial vector components as independent functions 
and obtain the eight other components JF, V and S^y as dependent functions. The 
results are: 

(a) transport equations: 

[A^ {n„^4] = {[A„p^ + n^],^^} - £,^^4A"A^ V'^] (D17) 
[A^ {n„ V^}] = {[A„p^ + n^], - £„^^4A"A^ A^] (D18) 

(b) mass shell equations: 

= [[P/. + ^^^.Pu + n.],in - [[A^, A,],i'^] 

+£a/3M4b" + n",A^],v^} (D19) 
M^v^ = [[p, + n^,p, + n,], - [[A,„ A,], vn 

+£a/3M.{[p" + n-,A^],i'^} (D20) 

(c) constraints: 

= [A^,V1 (D21) 



= {n^,>} (D22) 

= {n^, K} - {n,, V^} - £«;3^,[A", i^] (D23) 

(d) equations defining the dependent functions: 

2m:F^{U^,V^} (D24) 

2mV = [A^, >] (D25) 

2m4, = [A^, V,] - [A,, V^] + £a^^4n", A^}. (D26) 



APPENDIX E: VACUUM SOLUTIONS 
1. Massive quarks 

In tlie field- free limit the spinor equations can be solved exactly. Setting F^iy(x) = 
and choosing a gauge in which then also A^^x) = 0, the spinor equations reduce 
to 



4{E'-pl) + n 



4{E'-p',) + n 
4{E'-pl) + n 



V 



Q^-d^T-p^S^, 



(Ela) 
(Elb) 
(Elc) 
(Eld) 
(Ele) 
(Elf) 
(Elg) 
(Elh) 



2mV^ = 2p^T + d^S^, (Eli) 
2mA^ = -d^V + s^p^^p'S'^''. (Elj) 

Here E'^ = + m^, and the operator □ is defined through □ = d^d^. 

One can see easily that the transport equations (|El|i,ii,iii) are aheady contained 



in the constraints ([ET[vii,viii): For JF and V this follows simply by contraction with 
and respectively and exploiting the anti-symmetry of S^y = —Sy^. Eq. ([En|ii) 

follows from Eq. (^Ijviii) by contraction with e"' ^' p^' . 

Spinor components which vanish at t = — oo remain zero for all times: Formal 

integration of the transport equation Pfj,d^X{x,p) = for a spinor component X 

leads to 

PoX{t + At) = poX{t) - pid'X{t)At. (E2) 

With the initial condition X{t — oo) = we find X(t) = unless po = 0. But 
Pq = inserted into [4(i?^ — j9q) + 0]X = gives after formal integration 

doX{t + At) = doX{t) + (V^ - AE'^)AtX{t), (E3) 

so that with the initial condition X{t ^ — oo) = we have 

doX{t + At) =doX{t). (E4) 

Another formal integration then again yields X(t) = 0. Imposing for the spin and 
pseudoscalar densities vacuum initial conditions at t = — oo, <Sf^y{t — oo) = V{t 
— oo) = 0, we see that they remain zero for all times, and with them the axial vector 
density Afj, (see Eq. (pT|x)). 

The scalar density JF evolves according to 



d^P{x,p)={) (E5) 
[□ + 4(^2 - pI)]P{x, p) = (E6) 

lYioTn the first equation we get P{x,p) = P{p) while the second gives {E"^ — 
Po)^{p) = 0- The solution is stationary and homogeneous and can be written 
as a sum of particle and antiparticle contributions, 

Hx,P) = ^ [SiPo -E) + 5{po + E)] P{p), (E7) 

with an arbitrary momentum spectrum . The vector density is obtained through 
Eq. (^x) as 

V,{x,p) = ^:F{x,p) = [6{E~po) + 6{E+po)] Hp)- (E8) 



2. Massless quarks 

In the chiral limit the equations now read 



for al,^^ and 



for/W,SW 



= p^a^(±) (E9a) 
= a^a'^^^) (E9b) 
= d^af'^ - d,af^ ± 2£,^^,p"a^(±) (E9c) 



= 2p^/(±) - d^sfj (ElOa) 
= 9^/W + 2p'^SW (ElOb) 



With the help of the - hnearly dependent - "mass shell" equations we know that 
the solutions have to fulfill 

(V - n)^?^ = (Ell) 
(V - D)/^^^ = (E12) 
(V - n)^^) = 0. (E13) 

The solutions for the x dependence would therefore have the form exp{±p^x^^), which 
does not allow for reasonable boundary conditions at ±oo. Therefore we use 

the following Ansatz 

f^\x,p) = f^^\p)S{p') (E14) 
a^^\x,p)=a^^\p)5ip') (E15) 
S^^\x,p)=s^^\p)6ip'). (E16) 



Using Eqs. ( |E94|E10"a|) we find the additional restriction 



= ft^Sip,), /f ) = const. (E17) 
0'lf\x,p) =p^,ai^\p)5{p^) 

^(PO - \P\) - (^{PO + \P\) \ ^ P /-pisN 



Assuming now that the spin density vanishes for xq — oo we find that s^^^ has 
to vanish identically since s^^^ is independent of x^. The contribution from /*^^^ is 
restricted to = and constant in space and therefore irrelevant. Thus we are left 
with d^if \ giving contributions to vector and axial vector components as follows 



which have no longer a space and time dependence. 



APPENDIX F: MOMENT EQUATIONS IN SPINOR DECOMPOSITION 



The corresponding moment equations to Eqs. (4]^,|^) in spinor decomposition 
are for the constraints (Eq. (|4.3| )), which couple the n+V^ to all the lower moments, 



2 h+ilyo. 



(Fl) 



2 ["+ili'==-2m 



2 



2[n+l]^jfc. 



k=l 
_^Oijk 



(F3) 



2 ["+i]^o^_ 



2 [« + l]V,: 



2 ["+il4o=- 



:2m _ {iV(o) 0, NyO} - {n„ Wv^} _ 

-2m + ^"^^■'^[A,, ["^Vj-] - {N^Q^, '"U^'} + {^^ ["U°} + (F2) 

:2m - {iV(o)o , t"'^} - [A\ ["l5o.] - 

(2m ["U. + [A., '"^P]) - {iV(o)o, t"'^^''} 

{iV(o)o, t"U°} - {n„ ["U^} - 
EQ[%).,f"-'^'>]^+i 

-{iV(o)o, t"lV^} + {n„ ["IVo} + eo^,k[A\ ["U'=] - 
{iV(o)o, ["l5,o} - {n^ W^^.} + [A,, N^] + 



(F4) 

(F5) 
(F6) 



(F7) 



fc=i 



Oi\k 



(F8) 



and for the dynamical equations (Eq. ([4.4|)) containing explicit (9o derivatives 



k=l 



A;=l 



n 

E 

fc=i 

0=[A^,MV^] + E 



fc+i 



k=l 



o={n„ ["iv,} - {n,, Nv^} _ 5o.,,([A°, ["U'^] _ [A^ ["UO]) 



E 

A;=l 



fc+1 



fc+1 



-eo.,k ([M<i), [-'^U'^]^ - [Mf,), ["-'^U^ 
0=[Ao,["]^]-{n\W5o,} + 



E 

fc=i 



n 



fc+i 



0={n„ Wp} _ [A", N£„^] + [A^ M£^^.] + 



A;=l 



(F9) 



(FIO) 



(Fll) 



(F12) 
(F13) 



(F14) 



(F15) 
(F16) 



k+l 



Oi 



+ 



l\/Tj [n-k] ?• 



We used again the notations Eq. ( [4 .91) for the non-local operators M^^ and A^fc^^. 



APPENDIX G: THE COMMUTATOR EQ. ( [iTBl] ) 



We show that for classical external fields the commutator Eq. ( [4. 511 ) vanishes 
identically. Written out explicitly it reads 







2p^ + iD^{x) - ig J_^^^ ds{l - 2s) + tsdp)d; 

2p, - iD,{x) + ig r ds{l + 2s) + isdp)d'^ 

J-l/2 ^ 



P ' 



(Gl) 



To evaluate the commutator we first perform a Taylor expansion for the Schwinger 
string 



°° [isY 



X] 



n=0 ^• 
n=0 "'■ 

with the definitions Eq. ( |3.36| , p.3^ ). We will then need the two identities 



(G2) 



I3v 



n—1 



m=0 



n-l 



(G3) 



(G4) 



m=0 \m -r i-J 



The proof goes by induction as follows: 

The validity for n = is obvious. To make in Eq. (|G3|) the step from to n + 1 we 
need the relations 



V^AX -AV^X =[[D^,D,],X] = -tg[F^.,X] 



(G5) 



dp ■ ViAix), B{x)] = [{dp ■ V)A{x), B{x)] + [A{x), {dp ■ -D)B{x)] (G6) 



Thus we get 



Eq.jH) 



n-l 



m=0 



Eq. (M) 



n-l 



n 



m=0 



\{d,-vrF,^,{d,-vr-^Fp,]]d: 



-tg[F^a, (dp ■ VfFp.W + (9p • VY^^V^Fp, 



-^9 E L"! Am ■ vrF,^, id, ■ vr-^Fp^d; 



m=l 

n-l 



n 



m=0 V'''' -*-/ 

-igin + l)[F^^,id,-VrFp,]d; 

-ig[id, ■ WF,^, Fp,]d^ + idp ■ Vr+'V^Fp, 



m=l 



n 



,171, 



+ 



n 



, m + 1 , 



[id,-VrF,^,id,-V)"FpJ\ 



,m + 1 



n + 1 



id, ■ Vr+'V.F,, -zgY: L + 1 [(^P ■ (9, ■ Vr--F,J\d;. (G7) 



m=0 



This proves Eq. (|G3|) . For Eq. we find for n = 1 



[p^.d'^V^Fp,] = -g-^V^Ffs,. = -V^Fp,, 



(G8) 



i.e. tlie equation is correct. Tlie step n — ^ n + 1 is given by 



[p^, id, ■ Vf^'Fp,] 



[p^, id, ■ V)]id, ■ VTFp, + id, ■ V)[p„ id, ■ VTF,_ 



I3u\ 



Eqs. (|g4|,|g^) 



n 



n-l / ^ 



Eq.jH) 



m=0 



-{n + l){d,-'brV^Fp, 



n—l 



n 



n 



m=0 
n-l 

n-l / ^ 



\[{d,-vrF,^,{d,-vr-"^-'Fp^]d^ 



m=0 



{n + l){d,-VYV^Fp, 



m=l V"^/ 



m=l 
n-l 



m=0 



-{n + l){dp-VrV^Ffs, 



+ign[{d,-VT-^F^^,FpM 



n-l 



+ig Y: [{dp ■ Vr-'F,^, {d, ■ Vr-"^Ff,,]d; 



m=l 



^ ) + (m - 1)( ^ ) + m( ^ ^ 
ml \m I \m + 1. 



m 



n + 1 



+ IJ 

~{n + l){d,-'DY'D^Fp, 

+^9 E [{dp ■ Vr-'F,^, {d, ■ Vr--F,^]d^ □. (G9) 



m=0 



For the Schwinger string we get with the definitions of Eq. (|G2| ) the following 



commutator relations 



n=Q 



oo -n^n+l n 



n 



n n\ ' — 'n m + l\m 

n=0 m=U \ / 



[{d, ■ vrr,^, {d, ■ V)"FpJ\d;d^^ (GIO) 



and 



n=0 
oo -n „n+2 n 



n 



„=o ^0 + 2 



n + l 



In addition we find for tlie commutator of the Schwinger string 



(Gil) 



+ zsdp)d^, ^^^Fp,{x + ztdp)d^] = (G12) 

-E^ E {]s"^r-"^[{d,-vrF,^,id,-vr-"^F,^]dp^^. 



71 ' 

n=0 m=0 



Plugging in Eqs. ( plO| , |GTT| , |G12|) into the commutator Eq. (|G1|) and performing 
the s integration we get the following equation: 



-2ig{F^,{x) + F,^{x)) 



1 



n=0 
oo 



(n + l)(n + 3) 



V E E ■ i>rF,^. {dp ■ W-F.Jld^d^^ ^ 

n=0 m=0 \''V 



with 



9° 



(G13) 



f{n,m) 



1 



{n + 2)(n + 3) Vm + 1 
1 / 1 

+ 



+ 



2n + h 



n 



- m + 1 
2n + 7 



(n + 3)(n + 4) Vm + 2 n-m + 2 
2m + 3 

(m + l)(m + 2)(n — m + l)(n — m + 2) 



(G14) 



The first and second line on the r.h.s of Eq. ( |G13| ) vanish due to the antisymmetry 
of the field strength tensor and the Jacobi identity ( |3.18| ). In the third line one 
shows easily that for each combination of the indices n>0,0<m<n the function 
f{n, m) vanishes identically. Thus the commutator (|4.51|) vanishes. 

APPENDIX H: QED MOMENT EQUATIONS IN THE CLASSICAL LIMIT 

To derive the transport equations for the moments in leading order in the gradient 
expansion we need the moment equations to order h. These are with coupling 
constant g = —e inserted: 

2[n+i]vO(i)=2mW^(i) - 2piMv'(^) (HI) 

+^o.,7c ^^Q^ _ £5^F,,(x))Nvf + !^F,o(x)["-ilvf ) (H2) 

2[n+i]^(i)=2mWv°(i) - {d, - ^9^F,,(a;))W5°*(o) _ )[«-!] 50^(0) (^3) 

2[n+i]5iMi)=2m£°^^''=["Uf ^ _ 2p'^W50i(i) ^ 2p^T"l50^(i) (H4) 

2[n+i]^o(i)^_2p.M^*{i) (H5) 

2[n+l]yi(l)^2p*NV°^^^ 

+e''^' ((9, - ^9jF,Kx))N40) + ^F,o{xr-'W) (H6) 
2[n+i] 5»o(i)=_2pj. W 

+ (9, - £9^F,,(x))M^W + -F,o(x)["- (H7) 
c c 

2[n+i]p(i)=_(5. _ ^a^F,,(x))M/:0^(0) _ !^F,o(x)["-il/:°^(°) (H8) 

and 



-(9' - ^9p,F*^(x))MV°(°) - ^F,o(x)["-ilV°(°) (HIO) 

2^NV^(i)=2p,Wjr(i) + _ -diF,Jx)pS"'^''^ + — F.ofx)!"-^^''^^") (Hll) 

c ^ c 

2m["U°('^=2j9,["]/:°'(^) (H12) 

0=(9^ - -^di,F^,{x)pV^^''^ + ^F^o(x)["-^lV^(°) (H13) 

+5°*^'= - ^ajFfcKa;))["U°(°)^F,o(a;)["-iU°(°)) (H14) 

0=-2pi["]5°'(^) + (9o - -^d'pFoj{x)pJ^^^^ (H15) 

+ (9, - + — F,o(x)["-^lr^(°). (H16) 



To find the equations for tlie moments to zeroth order h we eliminate in Eqs. ( [Fl9| - 
[H16| ) witli n — > n + 1 tlie components of order h witli tlie lielp of Eqs. ([Hl[ - P^ ). 
Inserting Eqs. into Eq. (gT^) and Eqs. (g^) into Eq. ( [HT^ ) gives us (witli 



the unchanged Eq. ([H13|) ) the dynamical equations for the components Ny/^(o)^ 
c c 



+m(9* - + m^F'°(x)[""^l^(°) (H17) 



0=(9^ - -9^F^,(x))Mv^(°) + !^F^o(a;)["-i]V^(°) (H18) 



(H19) 



To apply now the constraints for the moments to zeroth order % 



m 



m 



[2n]y(0) ^ _g'2n P_[0] 

m 



[2n+l]y(0) ^ _g2n+l_P_[l] ^(0) 

Em 



(H20) 



[2n]^(0) ^^2n[0]j-(0) 



[2n+l] j^(0) ^ _g2n+l j-W 



we have to distinguish between even and odd n. For n even we find with 



pj 

'e 



(H21) 



from Eq. ( |H18| ) the following transport equation for the zeroth moment MjF^^^ 



= E' 



(H22) 



Similar we find for the odd moments with the help of i^^["lJ^(o) ^ Nj^(i) ^^j-^g transport 
equation for the first moment 1^1 JF^^^ 



= E' 



(H23) 



Expressed with fields E, B we finally get 



Q=E' 
0=E' 



(H24) 
(H25) 



The remaining dynamical equations Eq. ( |H17|jH19D prove to contain no new infor- 
mation: Eq. ([HT7D can be derived from Eq. ( [ET2^ ) or Eq. Eq. ( PTUD 
leads just to Eq. ( []i22| ) and Eq. ( |H23| ) with n replaced by n — 1. 

The transport equations for the other spinor components ^"'5^^-*, ^^'^If-* we find 
from Eqs. ( p9| , [Hll|JH14| , pl6|) for the n+ 1^* moments by eliminating the components 
to first order h with the help of Eqs. ( pl| - |H8|) . The transport equations are thus 



0=(9o - -c'^Fo,(x))["+iU°(°) 

c ^ c 

c ^ 



(H26) 



(H27) 



(5. - ^5ji^.Kx))["Um ^F,o(x)["-Uf 

_ -afcF*'=(x))["+iU°(°) - (^ + ^)V -°(x)["Uo(o) (H28) 
c c 

c 

c c 
-(9, - - i!i±ll^F,o(x)Wr^(o). (H29) 



Inserting the constraints for the moments to zeroth order h 



[2n]^0{0)^^2n^p . [0]^(0) 

m 



(H30) 



[2n]^(0)^^2nl ( [1]^(0) + ^ ^ (p 

m \ 



X 



[2n] s(^'>=~E^''p 



[i]cr(o) 



[2n]^{0)^^2n[0]^(0) 

1 [0]^(0) 
[2n+l]^0{0)^^2n+l_^p . 



[2n+l] ^{0)^^2n+l 



[0U(0) 



[2n+l] ^(0)^__gn2n+lp 



m \ E 



+ px [px 



E 



X 



E 



[2n+l]^{0)^_g2n+ 



[1U{0) 



(H31) 
(H32) 
(H33) 
(H34) 
(H35) 

(H36) 
(H37) 



we get for instance from Eq. (|H29|) for even and odd moments transport equations 
for the first two moments of (t^^\ Written with fields E and B these are 



0=E 



2n 



E{d' + -V,-E)-^ 



c 

e ( ( [iWW 



B X 



(H38) 



Q=E 



2n+l 



^(5° + -Vp ■ + (p . V + -(p X B) ■ Vp)— ^ 

c E 



c 

e / / Mo-W 



B X 



[1]^(0)' 



E 



(H39) 



The other dynamical equations do again contain no new information, leading to 
Eqs. ( [H38|JH39|) or px Eqs. ( [H38|jH39|) . 



REFERENCES 

[1] B. Miiller, The Physics of the Quark Gluon Plasma, Springer Verlag, Heidelberg, 
1985; 

J. Cleymans, R.V. Gavai, E. Suhonen, Phys. Rep. 130 (1986) 217. 

[2] Proceedings of Quark Matter 1993, Borlange, Sweden, Nucl. Phys. A 566 

(1994) , 

Proceedings of Quark Matter 1995, Monterey, California, Nucl. Phys. A 590 

(1995) , 

Proceedings of Quark Matter 1996, Heidelberg, Germany, Nucl. Phys. A 610 

(1996) . 

[3] M. LeBellac, Thermal Field Theory, Cambridge University Press, Cambridge, 
1996. 

[4] J. Kapusta, Finite Temperature Field Theory, Cambrigde University Press, 
Cambrigde, 1989. 

[5] E. Braaten, R. Pisarski, Nucl. Phys. B 337 (1990) 569; 

M. Thoma, in Quark Gluon Plasma 2, (R. Hwa, Ed.) World Scientific, Singa- 
pore, 1995, p. 51. 

[6] E. Laermann, Nucl. Phys. A 610 (1996) Ic. 

[7] Lattice 96, Proceedings, Nucl. Phys. B (Proc. Suppl.) 53 (1997). 

[8] S.R. de Groot, W.A. van Leeuwen and Ch.G. van Weert, Relativistic Kinetic 
Theory, North Holland, Amsterdam, 1980. 



[9] D.R. Nicholson, Introduction to Plasma Theory, Wiley, New York, 1983. 

[10] G.F. Bertsch, H. Kruse, S. Das Gupta, Phys. Rev. C 29 (1984) 673; 
H. Stocker, W. Greiner, Phys. Rep. 137 (1986) 277; 
G.F. Bertsch, S. Das Gupta, Phys. Rep. 160 (1988) 189; 
Gy. Wolf et. al., Nucl. Phys. A 517 (1990) 615; 

G. Batko, J. Randrup, T. Vetter, Nucl. Phys. A 536 (1992) 786; 

W. Gassing, V. Metag, U. Mosel, K. Nita, Phys. Rep. 188 (1990) 363. 

[11] P. Danielewicz, Ann. Phys. (N.Y.) 152, (1984) 239, and 305. 
[12] G.M. Welke, G.F. Bertsch, Phys. Rev. C 45 (1992) 1403. 
[13] S. Ochs, U. Heinz, Phys. Rev C 54 (1996) 3199. 
[14] U. Heinz, Phys. Rev. Lett. 51 (1983) 351. 

[15] U. Heinz, Ann. Phys. (N.Y.) 161 (1985) 48; and 168 (1986) 148. 

[16] E.P. Wigner, Phys. Rev. 40 (1932) 749; 

J.E. Moyal, Proc. Cambridge Philos. Soc. 45 (1949) 99; 

H. Grad, Commun. Pure Appl. Math. 2 (1949) 331. 

[17] P. Carruthers, F. Zachariasen, Rev. Mod. Phys. 55 (1983) 245 

[18] K.-C. Chou, Z.-B. Su, B.-L. Hao, L. Yu, Phys. Rep. 118 (1985) 235. 

[19] H.-Th. Elze, U. Heinz, Phys. Rep. 183 (1989) 81. 

[20] W. Botermans, R. Malfliet, Phys. Rep. 198 (1990) 115. 

[21] I. Bialynicki-Birula, Ann. Phys (N.Y.) 67 (1971) 252; 



M. Hillery, R.F. O'Connell, M.O. Scully, E.P. Wigner, Phys. Rep. 106 (1984) 
121. 

[22] H.-Th. Elze, M. Gyulassy, D. Vasak, Phys. Lett. B 177 (1986) 402; and Nucl. 
Phys. B 276 (1986) 706. 

[23] D. Vasak, M. Gyulassy, H.-Th. Elze, Ann. Phys. (N.Y.) 173 (1987) 462. 

[24] I. Bialynicki-Birula, P. Gornicki, J. Rafelski, Phys. Rev. D44 (1991) 1825 

[25] G. R. Shin, I. Bialynicki-Birula, J. Rafelski, Phys. Rev. A 46 (1992) 645; 
G. R. Shin, J. Rafelski, Phys. Rev. A 48 (1993) 1869; 
C. Best, J. M. Eisenberg Phys. Rev. D 47 (1993) 4639. 

[26] G. R. Shin, J. Rafelski, Ann. Phys. (N.Y.) 243 (1995) 65; 

C. Best, P. Gornicki, W. Greiner, Ann. Phys. (N.Y.) 225 (1993) 169. 

[27] P. Zhuang, U. Heinz, Phys. Rev. D 57 (1998) 6525. 

[28] F. Cooper, E. Mottola, Phys. Rev. D 36 (1987) 3114; and 40 (1989) 456. 

[29] J. Schwinger, Phys. Rev. 82 (1951) 664; 

V. S. Popov, Soviet Physics JEPT 35 (1972) 659; 

M. S. Marinov, V. S. Popov, Forts, d. Phys. 25 (1977) 373. 

[30] S. Ochs and U. Heinz in Thermal Field Theories and Their Applications, (Y.X. 
Gui, F.C. Khanna, Z.B. Su, Eds.) World Scientific, Singapore, 1996, p. 193. 

[31] M. Ploszajczak, M.J. Rhoades- Brown, Phys. Rev. Lett. 55 (1985) 147; 
M. Ploszajczak, M.J. Rhoadcs-Brown, Pliys. Rev. D 33 (1986) 3686; 
M.E. Carrington, M.J. Rhoades- Brown, M. Ploszajczak, Phys. Rev. D 35 (1987) 



3981. 

[32] P. A. Henning, Nucl. Phys. A 582 (1995) 633. 

[33] S. Mrowczynski, Ann. Phys. (N.Y.) 169 48 (1986); 

S. Mrowczynski, U. Heinz, Ann. Phys. (N.Y.) 229 (1994) 1; 
S. Mrowczynski, P. Danielewicz, NucL Phys. B 342 (1990) 345. 

[34] J. Schwinger, J. Math. Phys. 2 (1961) 407; 

P.M. Bakshi, K.T. Mahanthappa, J. Math. Phys. 4 (1963) 1, and 12; 
L.V. Keldysh, Sov. Phys. JETP 20 (1965) 1018. 

[35] R. Malfliet, in Correlations and Clustering Phenomena in Subatomic Physics, 
(M.N. Harakeh, J.H. Koch, O. Scholten, Eds.), NATO ASI Series, Vol B 359, 
p. 99, Plenum, New York, 1997. 

[36] H.-Th. Elze, Z. Phys. C 38 (1988) 211, and 47 (1990) 647. 

[37] H. Weigert, U. Heinz, Z. Phys. C 50 (1991) 195. 

[38] P. Zhuang, U. Heinz, Ann. Phys. (N.Y.) 245 (1996) 311; 
P. Zhuang, U. Heinz, Phys. Rev. D 53 (1996) 2096; 

A. Abada, M.C. Birse, P. Zhuang, U. Heinz, Phys. Rev. D 54 (1996) 4175. 

[39] Heribert Weigert, Diploma thesis. University of Regensburg, June 1989, unpub- 
hshed. 

[40] R. Fauser, Korrelationen in der Storungsreihe und in Transporttheorien fur 
Systeme im Nichtgleichgewicht, Ph.D. thesis, Ludwig-Maximilians-Universitat, 
Miinchen 1995, unpublished; 



P.A. Henning, R. Fauser, Boundary Conditions and Correlations in Path Inte- 
grals for Quantum Field Theory of Thermal and N on- Equilibrium States^ Los 
Alamos e-print archive |hep-pli/96(J5372 . 



[41] J. Hiifner, S. P. Klevansky, P. Zhuang, H. Voss, Ann. Phys. (N.Y.) 234 (1994) 
225; 

P. Zhuang, J. Hiifner, S. P. Klevansky, Nucl. Phys. A 576 (1994) 525; 

W. Florkowski, J. Hiifner, S. P. Klevansky, L. Neise, Ann. Phys. (N.Y.) 245 

(1996) 445; 

W. Florkowski, Description of hot compressed hadronic matter based on an ef- 
fective chiral lagrangian, Los Alamos e-print archive |hep-ph / 9 70 1 2 23 . 



[42] J.D. Jackson, Classical Electrodynamics, Second Edition, Wiley, New York, 
1975. 

[43] R. Hakim, H.Sivak, Ann. Phys. (N.Y.) 139 (1982) 230. 

[44] V. Bargmann, L. Michel, V. Telegdi, Phys. Rev. Lett. 2 (1959) 435. 

[45] S. Leupold, H. Weigert, Phys. Rev. D 54 (1996) 7695; 

S. Leupold, Feynman rules in radial gauge, submitted to Phys. Rev. D (Los 
Alamos e-print archive hep-th /9609222 ) . 



TABLES 







Am 


0*^ order 







1^* order 







TABLE I. The first two orders of the ?l-expansion for the non-local operators and 
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TABLE IL The first two orders for the gradient expansion for the non-local operators 
Va and for QED with external fields. 



